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APRESENTAÇÃO 


A coletânea “Estudios teórico-metodológicos en ciencias exactas, 
tecnológicas y de la tierra 2” é uma obra que visa identificar e solucionar 
problemas de qualquer área a partir de processos inovadores, e tem como 
objetivo principal a discussão científica por meio das obras que compõem seus 
capítulos. Este livro organiza trabalhos, pesquisas, estudos de campo e/ou 
revisões interdisciplinares, que integram e movem-se entre os diferentes ramos 
das engenharias, ciências e tecnologias inovadoras. 

O objetivo fundamental é classificar e apresentar com clareza as pesquisas 
teóricas-metodológicas realizadas em diversas instituições de ensino e pesquisa 
de diferentes países. A melhoria de processos organizacionais ou materiais, 
usando os avanços tecnológicos e modelagens matemáticas e/ou biológicas, 
levar-se-á a processos mais rápidos e eficazes em todos esses estudos. A 
necessidade da aplicação de novos conhecimentos relacionados à tecnologia 
em metodologias já conhecidas é importante, a fim de obter resultados práticos 
visíveis mais rapidamente, e esses resultados podem ser novos produtos, 
serviços ou processos, sempre desenvolvidos com o objetivo de alcançar algum 
tipo de avanço. 

Dessa forma, aqui são abordados diversos temas relevantes, com a 
finalidade de fundamentar o conhecimento da comunidade acadêmica e todos 
aqueles que, de alguma forma, se interessam pelas ciências exatas, tecnológicas 
e da terra em suas diversas áreas de aplicação. É crucial ter um material 
que demonstre progresso e aprimoramento em diversos processos, sejam 
eles organizacionais ou materiais, com informações que são extremamente 
relevantes, além de abordar temas atuais e relevantes para a sociedade. 

Assim, a obra Estudios teórico-metodológicos en ciencias exactas, 
tecnológicas y de la tierra 2 traz uma teoria sólida baseada nos resultados 
teóricos alcançados pelos vários professores, acadêmicos e pesquisadores 
que se dedicaram intensamente a desenvolver seus trabalhos, os quais serão 
apresentados de forma sucinta e didática aqui. Compreendemos a relevância 
da divulgação científica, portanto destacamos também a capacidade da Atena 
Editora em fornecer uma plataforma estabelecida e confiável para que esses 
pesquisadores possam apresentar e disseminar os resultados de suas pesquisas. 


Isabela Nascimento Souza 
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RESUMEN. La ingeniería de tejidos es un 
campo de investigación multidisciplinario 
que involucra biologia celular y molecular, 
ciencias de los materiales, bioingenieria, 
medicina, entre otras, que apuntan a 
restituir o reemplazar tejidos. Los enfoques 
predominantes incluyen el disefo y 
obtención de andamios que proporcionan el 
soporte estructural adecuado para el tejido 
óseo. En el presente trabajo, se propone una 
extensión del modelo matemático del tipo 
reacción-difusión propuesto por Svetlana V. 
Komarova et al. en 2003 de interacciones 
autocrinas y paracrinas entre osteoblastos 
y osteoclastos con términos biarmónicos 
que permite calcular la dinámica celular. Se 
realizó un estudio detallado de la dinámica 
y estabilidad del modelo de poblaciones 
estableciendo los valores de los parâmetros 
de control del modelo propuesto, bajo 
inestabilidad de Turing y Hopf. El modelo es 
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resuelto utilizando el método de Euler mejorado y diferencia finita en dos y tres dimensiones. 
Los resultados numéricos de la simulación computacional describen el proceso de 
remodelación ósea desde el punto de vista biológico. Posteriormente, mediante la técnica 
de proceso de disefio matemático y fabricación asistida por impresión 3D (MDP-3DPAM), 
se obtuvieron andamios que muestran una similitud con el tejido óseo trabecular humano, 
a base de ácido poliláctico (PLA), los cuales fueron caracterizados mediante pruebas de 
compresión mecánicas. Encontramos que la resistencia mecánica de nuestros andamios 
tiene un comportamiento viscoelástico que no se presenta en otros andamios obtenidos de 
la solución numérica del modelo de Komarova extendido por difusión empleando el mismo 
método. 

PALABRAS CLAVE: Remodelación ósea; Modelo matemático; Inestabilidad de Turing; 
Inestabilidad de Hopf; Andamios; Materiales biocompatibles. 


BIOLOGICAL MODELING FOR THE DESIGN, SIMULATION AND PROCESSING 
OF SCAFFOLDS FOR BONE REGENERATION BY 3D PRINTING 


ABSTRACT. Tissue engineering is a multidisciplinary research field that involves cell and 
molecular biology, material sciences, bioengineering, medicine, among others, that aim to 
restore or replace tissues. Predominant approaches include the design and procurement of 
scaffolds that provide adequate structural support for bone tissue. In the present work, an 
extension of the mathematical model of the reaction-diffusion type proposed by Svetlana V. 
Komarova et al. in 2003 of autocrine and paracrine interactions between osteoblasts and 
osteoclasts with biharmonic terms that allow calculating cell dynamics. A detailed study of the 
dynamics and stability of the population model was carried out, establishing the values of the 
control parameters of the proposed model, under Turing and Hopf instability. The model is 
solved using the improved Euler method and finite difference in two and three dimensions. The 
numerical results of the computational simulation describe the process of bone remodeling 
from the biological point of view. Subsequently, using the mathematical design process 
technique and 3D printing-assisted manufacturing (MDP-3DPAM), scaffolds were obtained 
that show a similarity with human trabecular bone tissue, based on polylactic acid (PLA), 
which were characterized by mechanical compression tests. We found that the mechanical 
resistance of our scaffolds has a viscoelastic behavior that does not occur in other scaffolds 
obtained from the numerical solution of the diffusion-extended Komarova model using the 
same method. 

KEYWORDS: Bone remodeling; Mathematical model; Turing instability; Hopf instability; 
scaffolding; Biocompatible materials. 


11 INTRODUCCIÓN 


En diversas ocasiones, distintas enfermedades de origen infeccioso, genético, 
degenerativo, agresiones físicas o químicas pueden dar lugar a una pérdida o alteración 
de las células de un tejido u órgano. Esta pérdida o dafio celular puede llevar a una 
alteración de la función normal del órgano y, como consecuencia, conducir al desarrollo de 


enfermedades o secuelas físicas que reducen la calidad de vida de una persona. Por este 
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motivo, uno de los principales retos de la medicina actual involucra la regeneración total 
y el restablecimiento de la función normal de los tejidos u órganos. Estos dos objetivos, 
regeneración y restablecimiento de la función normal de un tejido u órgano dafiado, son los 
fines principales de la medicina regenerativa e ingeniería de tejidos [1-2]. 

Actualmente la cirugía reconstructiva es el instrumento fundamental para la atención 
de estos pacientes; la utilización de órganos para trasplante usualmente se ve limitada 
por la escasez de donantes y el riesgo de transmisión de enfermedades. Esta necesidad 
de órganos y tejidos condujo a los investigadores a utilizar células para la reconstrucción 
de estos mismos. Hoy en día científicos de diversas áreas (molecular, celular, biológica) 
colaboran activamente con ingenieros mecánicos para desarrollar tejidos análogos que 
permitan a los médicos mejorar, mantener y restaurar la función de un órgano [3]. 

En el campo de la cirugía reconstructiva, se han generado soluciones basadas 
en trasplante de tejido xenogénico (origen animal), prótesis y manipulación localizada de 
sitios de lesión a nivel celular y molecular, en conjunto con los enfoques convencionales de 
reconstrucción a partir de tejido alogénico (origen en otro paciente humano) y auto génico 
(origen en el mismo paciente) [4]. 

La ingeniería de tejidos es un campo de investigación multidisciplinario que involucra 
biologia celular y molecular, ciencias de los materiales, bioingeniería, medicina, entre otras, 
que apuntan a mejorar o reemplazar funciones biológicas [5]. La ingeniería de tejidos 
se conecta con la práctica médica mediante terapias basada en células. Los enfoques 
predominantes incluyen métodos de recolección y expansión de tejidos vivos in vitro, así 
como el disefio de los biomateriales y andamios que proporcionan el soporte estructural 
adecuado para el tejido de interés. 

Uno de los objetivos más importantes de este campo es el disefio de los 
biomateriales que soporten las células, y de las moléculas de sefialización que permitan 
la regeneración del tejido, el cual lleve a una alternativa de trasplante de órganos y a la 
posterior reconstrucción de tejidos [6]. 

Uno de los principales problemas de los materiales utilizados para aplicaciones 
médicas es el déficit que presentan en cuanto a la analogía de propiedades (microestructura, 
composición, funcionalidad, biocompatibilidad, etc.) específicas del tejido a sustituir. En 
cuanto al perfeccionamiento de la estructura y propiedades de estos materiales, la biologia 
matemática puede servir como base fundamental, en donde tomando como base principios 
de bioquímica, física y fisicoquímica de los sistemas fuera de equilibrio [7-9] presentes 
en los sistemas biológicos y expresando mediante simulación; a partir de herramientas 
computacionales, la formación espacio-temporal y la dinâmica de estos procesos a partir de 
por ejemplo, el planteamiento de modelos deterministas de crecimiento [10] que permitan 
el observar, estudiar e incluso disefar la formación de tejidos u órganos a partir de la 
microestructura de la simulación resultante es muy prometedor. 


Un área interesante de estudio para la ingeniería de tejidos es el llamado proceso de 
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remodelación ósea. El remodelado óseo se produce en lugares espacial y temporalmente 
diferenciados del esqueleto e implica la reabsorción del tejido mineralizado existente 
por parte de los osteoclastos, seguida de la formación de hueso nuevo por parte de los 
osteoblastos [7, 11-13]. El modelado matemático de la remodelación ósea se ha centrado 
en varios aspectos de este proceso. Cada vez más, los datos experimentales sobre los 
mecanismos de sefialización entre las células óseas indican que las interacciones son 
complejas, y sacar conclusiones basadas únicamente en la intuición puede ser muy difícil y 
posiblemente engarioso. El modelado matemático proporciona una herramienta poderosa 
para reducir la ambigúedad en cuanto a causas y efectos en sistemas complejos, y permite 
probar varias hipótesis experimentales [14]. 

En la actualidad, los avances tecnológicos nos han permitido tener un crecimiento 
exponencial en el desarrollo de nuevos métodos de disefio y síntesis de materiales con 
potenciales aplicaciones en el área médica [3, 15, 16]. En cuanto a la síntesis de este 
tipo de materiales, en afios recientes con el avance de la tecnología de procesamiento de 
materiales, la impresión 3D resultó ser una herramienta bastante útil, la cual permite que un 
objeto tridimensional sea creado mediante la superposición de capas sucesivas de algún 
material deseado para un propósito en específico. 

Teniendo en cuenta este contexto, resulta hasta cierto punto sorprendente que la 
cantidad de modelos matemáticos que estudien los efectos integrados de las interacciones 
conocidas entre el hueso y las células propuestos a la actualidad sean pocos. Sin 
embargo, es cierto que ha habido algunos documentos recientes. Por ejemplo, los modelos 
matemáticos de la dinâmica de las células óseas que han sido propuestos por Komarova 
etal. [11] y Lemaire et al. [7]. 

En el presente trabajo, se propone un modelo matemático del tipo reacción-difusión 
(Komarova, 2003) de interacciones autocrinas y paracrinas entre osteoblastos y osteoclastos 
con términos biarmónicos que permite calcular la dinámica celular. Se realizó un estudio 
detallado de la dinámica y estabilidad del modelo de poblaciones estableciendo los valores 
de los parámetros de control del modelo propuesto, bajo inestabilidad de Turing y Hopf. 
El modelo es resuelto utilizando el método de Euler mejorado y diferencia finita en dos y 
tres dimensiones. Los resultados numéricos de la simulación computacional describen el 
proceso de remodelación ósea desde el punto de vista biológico. Posteriormente, mediante 
la técnica de proceso de disefio matemático y fabricación asistida por impresión 3D (MDP- 
S3DPAM), se obtuvieron andamios que muestran una similitud con el tejido óseo trabecular 
humano, a base de ácido poliláctico (PLA), los cuales fueron sometidos a pruebas de 
compresión mecánicas de donde encontramos que la resistencia mecánica de nuestros 
andamios es mayor al de otros andamios obtenidos empleando el método MDP-3DPAM. 
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21 MODELO EXTENDIDO KOMAROVA - SWIFT-HOHENBERG 


El modelo planteado por Svetlana V. Komarova et. al. en 2003 [11], es un modelo 
matemático que describe la dinâmica de la población de células óseas, con el número de 
osteoclastos y osteoblastos en una sola unidad multicelular básica (BMU) denotado por C 
y B, respectivamente. Las tasas de producción global de cada población celular reflejan 
el efecto neto del reclutamiento de precursores y la formación de células maduras. Las 
tasas de eliminación de células reflejan la muerte celular, así como la diferenciación de los 
osteoblastos en osteocitos y células de revestimiento del hueso. 

Se propuso que las células tengan la capacidad de interactuar entre sí a través de 
efectores, los cuales son liberados o activados por las células óseas y actúan de forma 
autocrina o paracrina de manera que, afectan localmente al tipo celular de origen o a otro 
tipo celular, respectivamente. En este modelo, tanto los osteoclastos como los osteoblastos 
pueden producir efectores locales capaces de activarse o inhibirse a sí mismos o a otro 
tipo celular. El ciclo de remodelación inicia cuando los osteoclastos se diferencian a partir 
de los preosteoclastos para obtener un aumento en la población, lo cual constituye una 
condición inicial para el sistema. Estos osteoclastos inician la reabsorción de hueso con 
la consecuente pérdida de masa ósea. Sin embargo, este evento a su vez provoca la 
diferenciación de osteoblastos a partir de los preosteoblastos, provocândose un incremento 
en la población de osteoblastos con la consecuente formación y restitución de masa ósea 
por acción de estas células. 

El modelo de remodelación ósea considera la región del hueso hasta el nivel celular, 
con influencia de factores locales y está basado en las ecuaciones: 


E = C9nB9n +(4c—B)C-Aç(1+V2)2C (1) 
+ = q9092B922 4 (Ap — B>)B — Ap(1+ v2)2B (2) 
= = —ky max(0,C(t) — €) + kymax(0,B(t) — B) (3) 


donde Cy B son el número de osteoclastos y osteoblastos; a, y B, son las actividades 
de producción y eliminación de células; los parâmetros g, representan la eficacia neta 
de los factores autocrinos o paracrinos derivados de los osteoclastos. Es importante 
mencionar que, una propiedad destacable es que las células son capaces de moverse 
de forma independiente a través de la matriz extracelular. Este movimiento celular puede 
definirse por la interacción de las células con sus vecinos más cercanos. Pueden sentir 
y responder más allá de las condiciones inmediatas del vecindario. Sin embargo, las 
células mesenquimales pueden responder a la interacción espacial de largo alcance de las 
células vecinas (efecto no local en la dispersión difusiva). Teórica y experimentalmente, la 
difusión puede ser reemplazada por la combinación de los mecanismos subyacentes con 
interacciones espaciales de largo y corto alcance, estos comportamientos biológicos se 
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consideran matemáticamente en los términos A, (1 + VC, A,(1 + V?)2B, con À, Actérminos 
constantes [22, 23]. 

Finalmente, Z representa la masa ósea de acuerdo con una tasa de recambio óseo, 
con las constantes k, y k, siendo la tasa de reabsorción y formación ósea respectivamente, 
por célula osteoprogenitora al día [11]. Sin embargo, esta suposición no es suficiente 
tomando en cuenta que nos encontramos en la descripción de un proceso con dinámica 
espacio-temporal. En cuanto a la resolución temporal del modelo descrito por Svetlana 
V. Komarova et. al. [11], toman en cuenta que el proceso de remodelación ósea es un 
fenómeno de dinámica temporal, el cual nos dará pauta para conocer las condiciones de 
comportamiento de evolución temporal del modelo. Los parámeiros k, y k, tendrían un 
significado de tal modo que la densidad de masa ósea disminuye o avanza con dinámica 
de acuerdo con los efectos netos de la reabsorción (C(t) > C ) y la formación (B(t) > B). Los 
parâmetros K,(t) y K,(t), como parâmetros dependientes del tiempo, satisfacen k (t) = rR, y 
ko(t) = r siendo R: 


= Jmax(o,c(t)-€) 


R= fmax(o,B(t)-B) (4) 


donde, es el período de los ciclos de C(t) y B(t), y r es determinada por la amplitud 
de las oscilaciones en cuanto a la evolución temporal de la densidad de masa ósea. El valor 
de R está bien definido, siempre y cuando B(t)B. 


81 PROCESO DE DISENO MATEMÁTICO Y FABRICACIÓN ASISTIDA POR 
IMPRESIÓN 3D (MDP-3DPAM) 

El estudio de las ecuaciones diferenciales se centra en la obtención de soluciones 
a partir de algunos métodos numéricos de ciertos tipos de ecuaciones y sistemas 
diferenciales, sin embargo, el obtener información cualitativa sobre los diferentes tipos 
de comportamientos o estados en las soluciones del sistema es vital para conocer el 
comportamiento del sistema de ecuaciones diferenciales. Es decir, averiguar qué tan 
estable o inestable es el sistema. A este análisis se le conoce como análisis de estabilidad 
lineal para sistemas de ecuaciones diferenciales parciales el cual es realizado en la etapa 
de modelado matemático y sirve como base para dar paso a la etapa del método numérico. 

El método numérico empleado para resolver el modelo extendido Komarova — Swift- 
Hohenberg, ecuaciones 1, 2 y 3, es el método pseudoespectral implícito descrito en [22] y 
esto se extiende al espacio 3D. El dominio espacial (discretizado por una malla de N, x Ny 
x Nz nodos) de esta ecuación es periódico en todas las direcciones (x, y, z). Este método 
basado en tales transformadas discretas brinda una precisión óptima en el sentido de que 
la tasa de convergencia para soluciones adecuadamente uniformes es más rápida que 
cualquier potencia de la malla de Nx x Ny x Nz nodos [21]. 

La metodologia para crear un archivo de disefio asistido por computadora (CAD) para 


Estudios teórico-metodológicos en ciencias exactas, tecnológicas y de la tierra 2 Capítulo 1 


un material polimérico con una distribución de tamafio de poro controlada se muestra en la 
Figura 1. La infografía muestra el proceso para desarrollar materiales con microestructura 
controlada. Partiendo de la solución numérica del patrón 3D obtenida del algoritmo del 
modelo extendido dado por las ecuaciones (1-3), el cual fue implementado en lenguaje de 
programación Python. La solución numérica en 3D es graficada la iso-superficie de mínimo 
valor en formato de archivo .pkl. Posteriormente se exporta para obtener un archivo con 
extensión .obj. Este archivo es redisefiado y puede traducirse a la extensión de lenguaje 
triangular estândar (.stl). En esta etapa, el tamafio de los poros se controla tanto como se 
desee. 


PROCESO DE DISENO MATEMÁTICO Y 


FABRICACIÓN ASISTIDA POR IMPRESIÓN 3D 


Etapa de Modelado Etapa de Disehio Etapa de 
Matemático Computacional Impresión 3D 


! 


Cc 22 
E CSM BI — (AÇO — Ac(1 + VC 


aB ape 
Gem meCDuBos — (, —As)B —Aa(1+ VB 


a =—ky max(0, C(t)- É) + ksmax(0,B(t) — 5) 


À Método numérico 


(FFT, pkl) Modelo 3 Rediseão (obj) ” 


Computacional Renderizado (.stl) 


Figura 1-. Diagrama de flujo de trabajo del proceso de disefio matemático y fabricación asistida por 
impresión 3D (MDP-3DPAM) para generar materiales porosos impresos en 3D con microestructura 
controlada. 


Los modelos 3D se pueden exportar a cualquier impresora generador de código 
de software para el paso final de procesamiento que también permite controlar el tamahio 
y los parâmetros de impresión. La extrusión de diferentes polímeros permite la creación 
de estructuras sólidas que poseen propiedades mecánicas complejas y una distribución 
controlada del tamafio de los poros. El proceso de impresión 3D que usa modelos 
matemáticos para disefiar el patrón [21] es un método completamente nuevo para disefiar 
y crear nuevos materiales avanzados. 


41 RESULTADOS Y DISCUSIÓN 


Soluciones numéricas del modelo matemático en 2D. 


Para los parámetros iniciales del modelo, se utilizó como estimación inicial 
lo ya propuesto por Komarova, que son datos experimentales obtenidos del análisis 
histomorfométrico de secciones óseas [11]. De esto se estimó que el início de la reabsorción 
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ósea en el hueso trabecular da lugar a la diferenciación de 10 a 20 osteoclastos que 
reabsorben hueso a un ritmo de 10 um/dia. 

En una misma localización, la erosión máxima se alcanza en un rango de 9 a 14 
días. Más tarde, unos 2000 osteoblastos Ilegan a la base de la cavidad y construyen hueso 
a un ritmo de 1 um/día, de manera que, de tres a cinco meses después de iniciarse el ciclo, 
el hueso vuelve a un estado quiescente (estado estacionario). 


— en-1 
Bi = 0.2 día Constantes de velocidad de 
B> = 0.02 dia? eliminación de células 


Valores iniciales de gij 


= cul 
q = 3 celulas día Constantes de velocidad de 
a, = 4 celulas dia”? formación de células óseas 


k, = 0.24% cell Idía”? Tasa de reabsorción y 
E formación ósea por célula 
k> = 0.0017% cell Idía”! osteoprogenitora al día 


Tabla 1. Resumen de parámetros usados en el modelo. 


Con estos parámetros y para la condición Z(0) = -k ()C + k(t)B (funciones) se Ilevó 
a cabo la simulación del modelo temporal de Komarova et. al. [11] utilizando el método 
numérico de Euler mejorado y diferencia finita en un arreglo cuadrado compuesto por una 
malla de 100x100 variando los parámetros de control g,, y À de acuerdo con los valores 
obtenidos del análisis de estabilidad lineal a 1,000,000 de iteraciones en un paso de tiempo 
dt=1x10. 

La solución numérica del modelo se obtuvo mediante la simulación computacional 
implementada en el software Phyton y sus librerías (pickle, Larray, pylab, Lgraphicarray, 
etc.), la programación del modelo fue elaborada en base al sistema de ecuaciones 
diferenciales parciales no lineales. Los resultados obtenidos de dicha solución numérica 
que muestran diferentes dinámicas de remodelación ósea presentes en el modelo con 
respecto a los cambios en la distribución de los osteoclastos y los osteoblastos se muestran 
a continuación. 

Caso 1: Donde 911 =—1.0, e =19y4; = 1.0. 
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Do O c) 


Figura 2-. a) Representa la densidad de biomasa ósea de la inestabilidad de tipo Turing, b) Representa 
la densidad poblacional de osteoclastos de la inestabilidad de tipo Turing, c) Representa la densidad 
poblacional de osteoblastos de la inestabilidad de tipo Turing. 


Caso 2: Donde 911 = 0.5, e = 1.8y4p = 1.0. 


Figura 3-. a) Representa la densidad de biomasa ósea de la inestabilidad de tipo Hopf, b) Representa 
la densidad poblacional de osteoclastos de la inestabilidad de tipo Hopf, c) Representa la densidad 
poblacional de osteoblastos de la inestabilidad de tipo Hopf. 


El requisito de estabilidad de un estado estacionario impone ciertas limitaciones 
de los valores que puede tomar cada parâmetro. Para la regulación autocrina de los 
osteoblastos, existe una solución estable con el comportamiento dinámico de un solo ciclo 
de remodelación cuando g,, es inferior a 0,9. Dado que el IGF es un conocido regulador 
autocrino positivo de los osteoblastos, para los cálculos se asume que g,, debe ser positivo 
pero pequefo. Un solo ciclo de remodelación es una solución estable para un amplio 
rango de valores positivos para la regulación paracrina derivada de los osteoclastos g,,, 
reflejando el hecho de que la activación de los osteoclastos conduce a un aumento de la 
formación de osteoblastos. La regulación paracrina g,, derivada de los osteoblastos se 
limita a valores negativos o positivos muy pequefios, lo que indica que el efecto neto de 
los factores producidos por los osteoblastos debe ser inhibitorio o neutro para formación de 
osteoclasto. En este sentido, los osteoblastos producen dos reguladores de osteoclastos 
opuestos, RANKL y OPG. El modelo sugiere que los efectos del RANKL no predominan. 

De los resultados obtenidos podemos ver que, sorprendentemente, los cambios en 
los factores locales autocrinos y paracrinos pueden determinar el tipo de remodelación 
ósea resultante. El factor autocrino g,, de los osteoclastos tiene la capacidad de cambiar el 
sistema entre modos de comportamiento dinámico. Una de las ventajas que tiene el haber 
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integrado los términos biarmonicos al modelo establecido por Komarova es que dichos 
parâmetros corresponden a un panorama biológico, mostrando una interacción celular, 
entre las células cercanas a su entorno y a largo alcance. 

En cuanto a los patrones que representan la densidad de biomasa ósea con la 
inestabilidad de tipo Turing y la inestabilidad de tipo Hopf (ver Figura 2 — 3) obtenidas 
tomando los parâmetros de control respectivos, muestra una distribución con máximos 
de densidad ósea [ver Figura 2 — 3 a), zonas blancas]; los cuales pueden ser atribuídos a 
zonas de hueso trabecular donde se encuentra el tejido óseo, y mínimos de densidad ósea 
[ver Figura 2 — 3 a), zonas negras] atribuidos a las zonas porosas presentes en el hueso 
trabecular (ausencia de tejido óseo). 

De estos resultados se pretende analizar los cambios de los patrones espacio- 
temporales que podrían surgir a partir de la modificación de parámetros de control claves 
en el proceso de remodelación ósea. Estos cambios tienen repercusiones biológicas en 
cuanto a la formación de diferentes tipologias e incluso patologías debidas a cuestiones 
bioquímicas y también dinâmicas del proceso. 


Solución numérica del modelo matemático en 3D. 


Respecto a la simulación 3D del modelo, nuevamente se utilizó el método numérico 
de discretización de laplaciano para la parte espacial y el método de Euler para la parte 
temporal mencionados anteriormente. Los valores de la tasa de reabsorción y formación 
ósea por densidad celular osteoprogenitora al día k (r, t) y k,(f, t) respectivamente con r=(x,, 
Y, Z,) para todo 1<1, j, k<100 para la solución de la densidad de biomasa ósea, son tomados 
como parámetros con evolución espacio-temporal para el cálculo de la microestructura 
tridimensional de densidad de tejido óseo. 

Las soluciones numéricas del patrón en 3D fueron obtenidas haciendo uso de la 
metodologia descrita en la sección 83. 


ce fa 


E. 115 


Mo. no 


o.105 

0.0995 
0.0944 
0.0893 
0.0843 


Figura 4-. a) Impresión en 3D del caso 1 a 1,000,000 de iteraciones en PLA (ácido poliláctico), bajo 
inestabilidad de tipo Turing, donde gu =—1.0, e 19y As =1.00, b) Solución numérica visualizada en 
Mayavi. 
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Figura 5-. a) Impresión en 3D del caso 2 a 1,000,000 de iteraciones en PLA (ácido poliláctico), bajo 


. a E Eu .z Bora E E 
inestabilidad de tipo Hopf, donde 9 =0.5, *=1.8 y As =1.0, b) Solución numérica visualizada en 


Mayavi. 


Analizando los resultados en 3D de dichos casos que representan tejido óseo sano 
de la inestabilidad de tipo Turing y la inestabilidad de tipo Hopf muestran una distribución 
con máximos normalizados de densidad ósea que van desde una zona donde existe una 
gran concentración de tejido óseo (100 por ciento) hasta zonas donde existe baja o nula 
concentración del tejido óseo. 


Procesamiento de andamios. 


Se procesaron 12 probetas con dimensiones de 25x25x25 mm tomando como base 
estructural las soluciones numéricas mostradas en las Figuras 4 — 5. Para cada caso fueron 
procesados 6 andamios, 3 de ellos bajo su posición original y los tres restantes con la 
posición rotada 90º hacia la derecha desde su posición original, esto con el fin de observar 
que la dirección de impresión no afecta a la estructura. Se utilizó una impresora 3D Select 
Mini V2 Monoprice para la obtención de los andamios, la temperatura de extrusión fue de 
210º a una velocidad de 50 mms. 

Para la realización de las pruebas mecánicas se sometieron las 12 probetas a 
ensayos de resistencia a la compresión con el equipo compresivo Universal Testing 
Machine CMSUTM-10 las cuales fueron determinadas mediante el software Smart Test, el 
equipo está configurado para soportar una fuerza máxima de 10 kN, con una pendiente de 
la rampa con valor de 3.00 mm/min, con condiciones de temperatura con valor de 25 grados 
Centígrados. Los parâmetros obtenidos del Módulo de Young, el límite de compresión y el 
límite de carga soportado por cada andamio se muestran en la Tabla 2. 
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Limite de == . 
ad cotagindai Límite de carga | Módulo de Young 
(MPa) (N) (MPa) 
Andamio 1.1-0 16.6 10302.71 35.5472 
Andamio 1.2-0 14.3 9081.51 32.18004 
Andamio 1.3-0 15.1 | 9287.61 35.2187 
Andamio 1.1-R 15.2 9328.31 29.1946 
Andamio 1.2-R 7.98 4966.0 23.9159 
Andamio 1.3-R 14.8 | 9225.51 33.1993 
Andamio 2.1-0 14.9 9251.41 36.1476 
Andamio 2.2-0 14.8 9251.81 39.9011 
Andamio 2.3-0 14.8 9270.81 32.2069 
Andamio 2.1-R 15.0 927911 28.6996 
Andamio 2.2-R 14.9 | 9215.11 39.1864 
Andamio 2.3-R 14.9 | 9266.41 36.6226 


Tabla 2-. Resultados de las pruebas de resistencia a la compresión de andamios basados en PLA. 


Las propiedades estructurales del hueso trabecular se determinan habitualmente 
mediante ensayos de compresión, tracción o flexión. Con los resultados obtenidos en estos 
distintos ensayos se ha observado que el hueso trabecular, de la misma manera que el 
cortical, presenta una mayor resistencia ante la carga de compresión que ante cualquier 
otro tipo [26]. La resistencia en ensayos de compresión tiende a variar entre 1,5 y 9,3 MPa 
en función de la región del esqueleto. 

Considerando los resultados podemos observar que el promedio de la resistencia 
en los ensayos de compresión de los andamios se encuentra en un intervalo de 7.9 a 
16.6 MPa. Por otro lado, el módulo de Young obtenido de las pruebas de resistencia a la 
compresión se encuentra en un rango desde 23 a 39 MPa. De acuerdo con lo reportado 
en la literatura, el módulo de Young para un hueso trabecular puede estar entre 10 a 500 
MPa, estas relaciones varían según los distintos tipos de hueso trabecular debido a las 
variaciones en la región anatómica, la edad y las enfermedades que alteran la arquitectura 
trabecular. 

Estos resultados se encuentran en un rango superior si los comparamos al 
establecido, sin embargo, este comportamiento puede ser atribuido a las propiedades 
mecánicas del PLA, lo cual es conveniente debido a que el PLA es un material que tiene 
propiedades de biodegradación en el cuerpo humano [29], brindándole mayor soporte 
durante el proceso de reabsorción y formación de hueso nuevo durante el proceso de 
remodelación ósea que sucede en el material. 


Los materiales termoplásticos o polímeros lineales suelen presentar una dualidad 


Estudios teórico-metodológicos en ciencias exactas, tecnológicas y de la tierra 2 Capítulo 1 


12 


de comportamiento ya que su respuesta de deformación frente a la aplicación de un estado 
tensional está constituida por dos componentes; una elástica la cual está relacionada con 
su naturaleza de sólido elástico y una viscosa, ligada a su naturaleza de líquido viscoso. 
Es por ello por lo que se considera que el comportamiento de estos materiales es de tipo 
“viscoelástico”, término que engloba ambas componentes en la respuesta del material [27]. 
Esta dualidad de comportamiento conlleva un cambio importante en el comportamiento 
mecánico, ya que sus propiedades mecánicas varían con el tiempo. 

Cuando un material plástico se somete a una tensión constante, se produce una 
elongación inicial inmediata, esto debido a la respuesta de sólido elástico y, además, 
se produce una elongación creciente relacionada con la componente de líquido viscoso 
del material. Cuanto mayor es el valor de la tensión aplicada, mayor es la elongación 
inicial elástica, y más rápido es el crecimiento de la elongación con el tiempo debido a 
la componente viscosa, este tipo de comportamiento se logra visualizar en las gráficas 
resultantes de los ensayos de compresión realizados a los andamios, comportamiento que 
se puede relacionar con la estructura del andamio obtenida de las soluciones numéricas 
del modelo matemático más que con el material, pues, realizando una comparación con 
los datos reportados por Fernández-Cervantes [18] (ver Figura 6), los cuales fueron 
reproducidos bajo las mismas condiciones que las soluciones numéricas propias, se 
logra observar una evidente mejora en la deformación que soportan los andamios al ser 
sometidos al ensayo de compresión, demostrando así que, el modelo matemático con una 
extensión de términos biarmónicos es eficiente en cuanto a la creación de andamios para 
la regeneración ósea. 


a) 1º b) | 
14 44 
” = 24 
E nd 
= 1 
e 19 E 4 
- e 1 
5 ê e sd 
ê E é] 
Ss: e 
a 44 ml Prueba 1 
—Pruebas 1 PruebaZ 
2 2 Prueba3 


oo 02 04 e 08 10 o 2 04 o) 08 10 
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Figura 6-. a) Gráfica de prueba de tensión-deformación, en color rojo, andamio basado en PLA con 
morfologia de tejido óseo del caso 2 con inestabilidad de tipo Hopf, en color negro, andamio basado 
en PLA con morfologia de tejido óseo de la solución numérica reportada por Fernández-Cervantes en 
[18] con inestabilidad de tipo Hopf, b) Gráfica tensión-deformación de andamios basados en PLA con 
morfologia de tejido óseo del caso 1 con inestabilidad de tipo Turing. 
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5 | CONCLUSIONES 


Se presenta un modelo matemático que examina las funciones cooperativas de la 
regulación autocrina y paracrina en el control de la remodelación ósea. El modelo matemático 
se basa en la suposición de que los efectores locales producidos por los osteoclastos y los 
osteoblastos regulan las tasas de formación de estos mismos, encontrando así que, el 
modelo predice diferentes modos de comportamiento los cuales se asemejan a la dinâmica 
de la remodelación ósea de tejido trabecular. 

El sistema es más sensible a la regulación autocrina de los osteoclastos, lo que 
refleja el hecho de que los osteoclastos reabsorben el hueso como pequefios equipos de 
células muy activas, que se reclutan rápidamente y luego se eliminan. Contrario a eso, los 
osteoblastos son mucho menos activos y los cambios en el número de osteoblastos se 
producen con mayor lentitud; en consecuencia, se necesitan muchos más osteoblastos en 
un mismo lugar de remodelación ósea. 

Del procesamiento de los andamios se comprobó que la técnica de proceso de disefio 
matemático y fabricación asistida por impresión 3D (MDP-3DPAM) permite la reproducción 
exitosa de piezas con morfologia de tejido óseo. Los andamios obtenidos son comparados 
con el tejido óseo trabecular humano sano, en donde en ambos casos existen trabéculas 
interconectadas. Los resultados de las propiedades mecánicas mostraron una similitud 
con respecto a las propiedades mecánicas del hueso trabecular humano, así mismo, se 
observa que cuando los andamios están sometidos a una deformación constante con el 


paso del tiempo, estos dan lugar a un fenómeno ligado a viscoelasticidad. 
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CAPÍTULO 2 


LOS TEOREMAS DE LA FACTORIZACIÓN DE 
CORDERO EN EL CONJUNTO DE LOS NÚMEROS 
ENTEROS 


Ronald Cordero Méndez 
Universidad Internacional San Isidro 
Labrador. Costa Rica. 


RESUMEN: Se presentan Los Teoremas 
de la Factorización de Cordero en el 
conjunto Z. Ejemplos de aplicación de 
estos teoremas. Utilidad en el cálculo de 
números primos muy grandes, así como 
en la factorización de números enteros 
muy grandes. Contribución a la Teoría de 
Números en un problema del que no se 
tiene una solución eficiente como es la 
factorización de números enteros y material 
de apoyo para la construcción de software 
que permitan el cálculo de números primos 
muy grandes necesarios en la criptografia. 
PALABRAS CLAVE: Factorización, 
números primos, números afortunados de 
Euler, Teoremas. 


ABSTRACT: The Cordero Factorization 
Theorems in the set Z are presented. 
Examples of application of these theorems. 
Useful in the calculation of very large prime 
numbers, as well as in the factorization of 
very large integers. Contribution to the 
Theory of Numbers in a problem for which 
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there is no efficient solution, such as the 
factorization of integers and support material 
for the construction of software that allows 
the calculation of very large prime numbers 
necessary in cryptography. 

KEYWORDS: Factorization, prime 
numbers, Euler's lucky numbers, Theorems. 


11 POLINOMIOS GENERADORES 
DE NUMEROS PRIMOS Y 
COMPUESTOS. 

Los números primos han sido 
tema de muchas investigaciones, muchas 
repetitivas que contribuyen poco al tema, 
lo que verifica la frase del gran matemático 
Leonhard Euler que dice: “Los matemáticos 
han intentado en vano, hasta la actualidad 
descubrir algún orden en la secuencia de 
números primos, y tenemos razones para 
creer que se trata de un misterio que la 
mente humana nunca resolverá” (Leonard 
1707-1788, 
Camacho y Camacho, 2020, p.85). Hasta 


Euler, mencionado | por 
el momento en el afio 2022 este misterio 
no ha sido resuelto, por lo que creo que 
Euler puede estar en lo cierto. Leonhard 
Euler nació el 15 de abril de 1707 en 
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Basilea, Suiza y murió el 18 de septiembre de 1782 en San Petersburgo, Rusia (Aznar, 
2007). Extraordinario matemático del siglo XVIII. 
Otra frase que lo afirma dice: 


El encanto de los números primos consiste quizás en la imposibilidad 
de explicar en qué orden aparecen. Cada uno se dispersa a su antojo, 
cumpliendo la condición de no tener más divisores que el uno y él mismo. 
Aunque no cabe duda de que cuanto más grandes son, más difícil resulta 
encontrarlos, y es imposible predecir su aparición siguiendo ninguna regla... 
"La fórmula preferida del profesor (Ogawa, 2003, mencionadado por Frases y 
Pensamientos, s.f., párr. 4) 

Nuestra pregunta ahora es cómo encontrar números primos, si no es posible 
encontrar una fórmula polinomial o de otro tipo que nos genere todos y cada uno de los 
números primos, o por lo menos una fórmula que genere solamente números primos 
aunque no sean consecutivos. 

En algún momento dado aparecen los números compuestos que se mezclan con los 
números primos, por lo que me lleva a suponer que el cribado es una buena opción para 
encontrar números primos grandes. 

Con ayuda de los polinomios P(n) =n? + n+ p, donde p=2,3,5, 11,17, 41 que 
resulta ser polinomios que generan números primos cuando n toma valores desde O hasta n 
= p-2,y luego generan números compuestos y primos mezclados, por lo que el problema 
de encontrar una fórmula que genere solamente números primos no lo resuelve este tipo 
de polinomios. Pero encontrar una fórmula que genere los números compuestos que 
son generados por estos polinomios es el tema de la investigación además de buscar un 
procedimiento que ayude a cribar los números primos. 


21 POLINOMIOS DE LA FORMA P(n) = nº + n+ p, DONDE p=2,83,5, 11, 17,41. 


Los polinomios P(n) = nº + n+ p, generan números primos, por ejemplo, se generan 
los números primos: 41, 43, 47, 53, 61, 71, 83, 97, 113, 181, 151, 173, 197, 223, 251, 281, 
313, 347, 383, 421, 461, 5083, 547, 593, 641, 691, 743, 707, 853, 911, 971, 1033, 1097, 
1163, 1231, 1301, 1373, 1447, 1523,1601 

cuando P(n) = nº? +n+41y desde n=0Ohastan=41 -2=39, en total 40 números 
primos, pero a partir de n = 40 se generan números compuestos y números primos. A este 
poliomio se le Ilama polinomio de Euler. 

Otra forma de escribir el polinomio de Euler es P(n) = nº - n+ 41, pero éste genera 
los números primos anteriores cuando, n toma valores desde 1 hasta 40. 

Otro polinomio de esta forma que genera números primos es P(n) ="? + n+ 17 
desde n= 0 hasta n= 17 — 2 = 15, el cual fue descubierto por el matemático Adrien Marie 
Legendre: 
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Legendre nació en París en el afo 1752 en una familia rica. Recibió educación 
en el Collage Mazarin en París, y defendió su tesis en física y matemática 
en 1770. Murió en París en el afão 1833, después de una larga y penosa 
enfermedad. Su viuda conservó cuidadosamente las pertenencias del 
matemático para preservar su memoria. El último lugar donde vivió fue en el 
pueblo de Auteuil en París, Francia (Fernández y Tamaro, 2004, párr.1) 


31 LOS NÚMEROS AFORTUNADOS DE EULER 


Primero dejemos claro que Goldbach y Legendre demostraron que no es posible 
encontrar un polinomio que dé números primos para todo número natural, el primero lo 
demostró para coeficientes enteros y el segundo para funciones algebraicas racionales. 

El matemático Rabinowitz demostró que P(n) = nº + n + p da números primos para 
n=0,...,p-2 siy solo si 1 - 4p es el negativo de un número de Heegner, que son los 
únicos números positivos k, que cumplen no ser cuadrados perfectos y que en el anillo de 
enteros del cuerpo Q(v-k) es de factorización única. 

Los números de Heegner son: 1,2,3,7,11,19,43, 67,168. 

Además los números afortunados de Euler son los enteros positivos p, para los que 
1 -4p =—k, siendo k un número de Heegner, y mediante comprobación obtenemos que los 
únicos posibles son 2, 3, 5, 11, 17, 41 y el número de Heegner asociado al 41 es el 163. 


41 PROBLEMÁTICA DE LA FACTORIZACIÓN DE LOS NÚMEROS ENTEROS. 


No se ha encontrado un procedimiento eficiente para factorizar completamente (en 
sus factores primos) un número entero muy grande. Los números enteros muy grandes 
más difíciles de factorizar son aquellos que tienen solamente dos factores diferentes de uno 
y de aproximadamente el mismo tamafio. 

En 1970 el límite era 20 dígitos, luego con el método de las fracciones continuas de 
Brillhart- Morrison se Ilegó a 50 dígitos en 1980. En el ao 1994 se logra factorizar números 
de 100 dígitos por medio del método de la criba cuadrática de Carl Pomerance y luego se 
afiaden métodos de factorización de números enteros como el de cuerpos algebraicos de 
Pollard y el método de Lenstra el cual utiliza curvas elípticas. 

El algoritmo de criba cuadrática (QS), es un algoritmo de factorización de números 
enteros y en la práctica es el segundo método más rápido conocido superado en poco por 
la criba general del cuerpo de números. Es de propósito general, lo que significa que su 
tiempo de ejecución únicamente depende del tamafio del número entero y no sobre una 
estructura especial o de propiedades. 

La factorización de curva elíptica de Lenstra o método de factorización de curva 
elíptica es un rápido algoritmo de tiempo de ejecución sub-exponencial para la factorización 
de enteros que utiliza curvas elípticas. Es el tercer método más rápido conocido de 
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factorización, superado por el algoritmo de la criba cuadrática (QS) y por la criba general 
del cuerpo de números. 

El método de factorización de curva elíptica es de propósito especial y el más 
adecuado para encontrar factores pequefios, útil para encontrar divisores que no superen 
los 20 o 25 dígitos, así como su tiempo de ejecución está dominado por el tamafio del factor 
más pequefo, en lugar de por el tamafio del número a factorizar. Frecuentemente este 
método se utiliza para eliminar factores pequefios de un entero muy grande con muchos 
factores; si el entero resultante es todavía compuesto, entonces solo tiene factores grandes 
y es factorizado mediante el uso de técnicas de propósito general. El factor más grande 
encontrado con este método cuenta con 75 dígitos y fue descubierto el 2 de agosto del 
2012 por Samuel Wagstaff. Incrementando el número de curvas probadas se mejoran las 
posibilidades de encontrar un factor. Hay consenso de que la factorización de números es 
un problema difícil, pero los avances nos dan una esperanza, quizás haya una solución 
muy pronto. 

Hoy día la computación cuántica enfrenta el reto matemático de factorizar en números 
primos números enteros muy grandes. Investigadores de la Universidad Politécnica de 
Madrid en Espafia han abordado el reto de factorizar números muy grandes en sus factores 
primos utilizando un dispositivo cuántico que simula la aritmética, en vez de calcular. 

Pero la complejidad para factorizar un número muy grande se ha aprovechado para 
realizar algoritmos de criptografia que permiten mantener información privada segura. La 
criptografia es el estudio de las técnicas matemáticas relacionadas con los aspectos de 
seguridad informática tal como: la confidencialidad, la integridad de datos, la autenticidad 
y el no rechazo. Hay dos tipos de criptografia: criptografia simétrica o de clave privada 
y criptografia asimétrica o de clave pública. La principal aplicación de la criptografia es 
proteger la información para evitar que sea accesible a observadores no autorizados. 
Verificar que un mensaje no haya sido modificado a conveniencia por ejemplo. Es aquí 
donde la matemática juega un papel muy importante, y sobre todo los números primos 
muy grandes, puesto que los números primos muy grandes de más de 100 dígitos permiten 
formar semiprimos o biprimos que son el producto de dos números primos muy grandes y 
que sirven de códigos de seguridad, por su dificutad para ser factorizados. 

La factorización de un número natural o entero necesita de más estudio, actualmente 
no se tiene un procedimiento 100% eficaz y eficiente, y esto no es producto de falta de 
investigación, puesto que en la historia de la matemática podemos considerar que valiosos 
genios de la humanidad trabajaron para encontrar una solución, como es el caso del francés 
Pierre de Fermat (1601-1665) y el Suizo Leonhard Euler (1707-1783). 

El creciente desarrollo de la tecnologia y la construcción de diferentes algoritmos 
matemáticos no han podido darle fin a la rapidez y eficiencia que debe tener un software 
que factorice números enteros, esto deja en evidencia el reto de los matemáticos actuales 


y futuros en el estudio de este tema. 
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51 NÚMEROS COMPUESTOS Y PRIMOS GENERADOS POR POLINOMIOS. 


é Puede un polinomio no constante, de coeficientes enteros, tomar solamente valores 
primos”; la respuesta es no, é; por qué?, debido al siguiente teorema. 

Teorema. 

Si f(x) € Z[x] tiene grado positivo, entonces existe una cantidad infinita 

de números naturales n para los cuales f(n) es un número compuesto. 

La otra pregunta que podemos hacernos es la siguiente: 

é Todo polinomio no constante f (x) E Z[x] (conjunto de todos los polinomios de 
coeficientes enteros) puede generar infinita cantidad de números primos? Intrigante 
pregunta. Bouniakowsky y luego Schinzel y Sierpinski (1958) conjeturaron que si f (x) E 
ZIx] irreducible, primitivo, es decir, el máximo común divisor de sus coeficientes enteros 
es igual a 1 y si más aún no hay ningún primo que divida a todos los valores de f(n) para 
enteros arbitrarios n, entonces este tipo de polinomios siempre generan números primos. 
(Revista Colombiana de Matemáticas. Vol XXI (1987). pág 263). 


61 EL TEOREMA PRINCIPAL. 


TEOREMA. 
Sea p un número primo y hagamos E =m"?+n+ op. Entonces las siguientes 
condiciones son equivalentes: 
IE p=2,3,4,11,17,41 
ii. f, (n) es un número primo para n=0,1,2,3,...,p-2 


ii. Q(/1 — 4p) tiene número de clase igual a 1. 


Tenemos que i. — ii., basta con sustituir p por 2,3,4,11,17,41 y verificar que en efecto 
se generan números primos desde n=0 hasta n=p-2. 

La equivalencia de las condiciones ii. > iii. fue demostrada por primera vez por 
Rabinovitch en 1912 y en 1936 es demostrada por Lehmer. Mientras que iii. > ii. es 
demostrada por Szekeres en 1974 y por Ayoub y Chowla en 1981. (Revista Colombiana de 
Matemáticas. Vol XXI (1987). pág 281). 


7 | IMPORTANCIA DE FACTORIZAR NÚMEROS ENTEROS GRANDES 


Los algoritmos criptográficos se basan en un un problema computacionalmente difícil 
de resolver, la mayoría de los algoritmos empleados en criptografia de Ilave pública, tienen 
como problema a resolver la factorización prima de números grandes. El más usado en la 
actualidad es el algoritmo RSA, descubierto por Ron Rivest, Adi Shamir y Len Adleman, de 
ahí sus iniciales RSA. El RSA es el algoritmo de llave pública mas usado en el planeta, 
tanto para cifrado y descifrado como para fines digitales. 

El algoritmo RSA es seguro mientras la Ilave privada solo la conozca el duefio. La 
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llave pública y privada se encuentran matemáticamente relacionadas de tal manera que a 
partir de la llave pública sea imposible generar la llave privada. 
La fortaleza del algoritmo radica en la dificultad para factorizar números grandes, y 


lo cual hasta el momento no tiene una solución satisfactoria y definitiva. 


81 LOS TEOREMAS DE LA FACTORIZACIÓN DE CORDERO EN EL CONJUNTO 
DE LOS NUMEROS ENTEROS. 

Los teoremas aquí publicados, y sus aplicaciones, aunque no resuelven el difícil 
problema de factorizar completamente un número entero muy grande, son muy útiles para 
factorizar en dos factores un número entero, sin necesidad de recurrir a software, test de 
primalidad, o divisiones sucesivas. La factorización de estos números de forma polinomial 
cuadrática se puede encontrar utilizando fórmulas matemáticas descubiertas por el autor 
de este artículo. 


8.1 Primer Teorema de la factorización de Cordero en Z 


Teorema que permite factorizar en dos factores números de la forma ?P + n+p 
donde ne Z 
Sear,s,x,k EZ y pe(2,3,5,11,17,41). 


D Si n=((rs-r+12ex2?-Gs-r+(rs-r+2s-Dex+tpe(rs —-r+ 1)2+ 
(rs+s+HD(s-D+r+2)Kk-D+r(rs-r+D)ex?—- (r(r+2)(s—-2D) +72 +37 + 
Dx+pr(rs-r+D+(r+1)(s—1) 


entonces: n2+n+p es compuesto y dos de sus factores tienen la forma: 


fh=(rs-r+1)2ex2-(rs-r+D(rs-r+2s—-Dex+pe(rs—-r+ D+(rs+s 
+D(S-D+r+2) 


n2+n+p 


f= E A*K-DAMLADK- Der -r(r+2x+pr+r+1 


2) Si n=((rs-r+1)ex2-(rs-r+D(s-r+2s-Dex+pe(rs —-r + 1)2+ 
(s+s+HD(s-D+r+ZD)KAD—(r(rs-r+1)ex2?—- (r(r+2)s— 2) +72 +3r+ 
Dx+prírs-r+D+(r+U(s-D)+1) 


entonces: n?+n+p es compuesto y dos de sus factores tienen la forma: 


h=(s-r+12ex-(rs-r+1)rs-r+2s-Dex+pe(rs-r+1)+(rs+s 
+1)(s—- 2) +7r+2) 


n2en+p 


fr = A 


Con B=r(rs-r+1D)ex2?-(r(r+4+2)s-2D+7r2+3r+Dx+tprírs-r+D+(r+D(s—1) 


fKAD CADA D+ArS -r(r+2Dx+pr+r+1 


Nota: 17 y sno pueden ser cero simultáneamente. 


También es importante aclarar que el teorema es válido en el conjunto de los 
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números reales y en el conjunto de los números complejos, es decir que las variables 


pueden ser sustituidos por números complejos y el Teorema se cumple. 
8.1.1  Aplicaciones del Primer Teorema de La Factorización de Cordero 
en el Conjunto Z. 


Si damos valores enteros a las variables que conforman las fórmulas del Primer 
Teorema de La Factorización de Cordero en los enteros, podemos encontrar la factorización 
en dos factores de los números de la forma n2 + n+ pconp € (2,3,5, 11, 17, 41). Es 
importante aclarar que los cinco teoremas que aquí se publican son válidos en cualquier 
conjunto numérico, es decir las variables pueden ser sustituidas por números naturales, 
enteros, racionales, irracionales, reales o complejos y los teoremas se cumplen. Pero el 
objetivo de esta publicación es factorizar números naturales, por lo que nos limitaremos a 


dar solamente valores enteros a las variables que contienen las fórmulas. 
Aplicación 1 
Sea r=239,s=678,x=100,k=451yp=t41 tenemos que: 
rs—r+1=239+678-239+1 = 161804 
rs —- r+2s-1=239+678-2394+2+678-1= 163158 
rs+s + 1=239+6784678+1= 162721 


Luego: 


h=(s-r+19).x22-(s-r+Das-r+25-Dex+tpe(rs-r+1)+(s+s 
+D(s-D+r+2) 


f, = 1618042 x 1002 — 161804 x 163158 * 100 + 41 * 1618042 + 162721 + 676 
+ 241 = 260238894367493 


B=r(rs-r+1)ex2-(r(r+2)(s-2) +72 +3r+1)x+pr(rs-r+1) 


+(Gr+1)(s-1) 


B = 239 + 161804 + 1002 — (239+241+676+ 2392 +3+239+1)+100+41+239 
* 161804 + 240 + 677 = 384397763576 


n=fh+(k-D+B 


n = 260238894367493 + 450 + 384397763576 = 117107886863135426 
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b=f,*(K-1)2+(Q8+D(k-D)+rixor(r+2)x+pr +r+1 


fa = 260238894367493 + 4502 + (2 * 384397763576 + 1) * 450 + 2392 + 1002 — 239 x 241 
+ 100 + 41+2392 + 239 + 1 = 52698722067972343651 


Por el Primer Teorema de la Factorización de Cordero en los enteros, 


117107886863135426? + 117107886863135426 + 41 
= 260238894367493 + 52698722067972343651 


Factorizado en sus factores primos es: 


117107886863135426? + 117107886863135426 + 41 
= 4943 * 150401 + 16419517 + 649 305521 * 1730 300293 


Aunque el Teorema no logra la factorización completa, es muy útil porque logra 
factorizar el número en dos factores utilizando fórmulas matemáticas y sin necesidad de 
recurrir a test de primalidad, calculadoras factorizadoras de números, divisiones sucesivas 
o software. 

Obsérvese que a pesar de que: 

1171078868631354262 + 117107886863135426 + 41 = 137142571655489270032 
49522479336943 es un número grande (36 dígitos) su factorización prima tiene solamente 
5 factores, 


Aplicación 2. 
Sea r = 139, s=1608x=-15,k=-451 yp= 11 tenemosque: 


rs — r + 1=139+1608-139+1 =223374 
rs-— r + 2s-1=139+1608-139+2+1608-1= 226588 
rs+s + 1=139+1608+ 160841 = 225121 
Luego: 


h=(s-r+19ex-gs-r+Ds-r+2s-Dex+pe(rs-r+1)+(s+s 
+1)(s-2)+r+2) 


fi = 2233742 + (—15)? — 223374 + 226588 * (—15) + 11 + 2233742 + 225121 
* 1606 + 141 = 12535012317883 


B =r(rs-r+1).x-(r(r+2)(s-2) +72 +3r+1)x+pr(rs-r+1) 


+(r+1)(s-1) 


B = 139+223374+ (15)? — (139 + 141 * 1606 + 1392 +3+139+1)+*(-15) + 11+139 
+ 223374 + 140 + 1607 = 7800221671 


n=f*(k+D-(B+1) 
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n = 12535012317883 + (—-450) — 7800221672 = —5640763343269022 


fR=ftKAD- MADURA D+Hrgr(r+2Dr+tpr+r+1 


f, = 12535012317883 + (—450)2 — (2 * 7800221671 + 1) * (—450) + 1392 x (—15)2 — 
139+141+(—15)+11+1392+139+1 = 2538347014575665731 


Por el Primer Teorema de la Factorización de Cordero en los enteros, 


(-5640763343269022)2 + —5640763343269022 + 11 
= 12535012317883 + 2538347014575665731 


Factorizado en sus factores primos es: 


(-5640763343269022)? + —5640763343269022 + 11 
= 1122 + 797 * 129981359 + 245 154001 * 10354 091731 

Su factorización completa tiene 6 factores que son pocos, tomando en cuenta que el 
número a factorizar es grande. Utilizando calculadoras factorizadoras de números enteros, 
como la de Dario Alpern (véase, 4. Factorization using the Elliptic Curve Method) de https:// 
www.alpertron.com.ar/ECM.HTM. Podemos factorizar números muy grandes con factores 
que sobrepasan los 100 dígitos, véase www.REVISTAELLABRADOR. Primer Teorema de 
La Factorización de Cordero en los números enteros. 


8.2 Segundo Teorema de la factorización de Cordero en Z 

Este teorema es la culminación de un resultado matemático (El teorema de la 
multiplicación de Cordero en Z), publicado por Ronald Cordero Méndez en: Memorias 
de Fimat 2020, que lo puede encontrar en: libro memorias fimat concites 2020 c.pdf, 


página 104. 
Sea sx, TEZ,s*0, pe(2,3,5,11,1741) 


1) Si n=(sx2+s(s—-2)x+ps?-s+1)*T+sx2+(s— 1)x+ps—1 entonces 
n?+n+p es compuesto y dos de sus factores tienen la forma: 
fi=six+s(s—-2)x+ps?-s+1 


n2+n+p 2 2 2 
ja = = +28 +1)T+x+x+p conB=sx+(s-I)x+sp—1 
1 


2) Si n=(stx2+s(s—-2)x+ps?-s+1)+*T-—(sx2+(s—1)x+ps) entonces n? + 
n+p es compuesto y dos de sus factores tienen la forma: 


fi=sêx+s(s-2)x+ps?-s+1 


h= COM = AT? (2B+NT+R+z+p conB=sx+(s-1I)x+sp-—1 
1 
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Este Teorema al igual que el primer teorema, su objetivo es encontrar una 
factorización en dos factores de números de la forma n? +n + p con p € (2, 3,5,11,17,41). 


8.2.1 Aplicaciones del Segundo Teorema de la Factorización de Cordero 
en los números enteros. 


Aplicación 1 
Sea s=678,x= -—-80, T=12 yp = 41 tenemos que: 


fh=sx+s(s—-2)x+ps?-s+1 
fi = 6782(-80)? + 678 +676+ (-80) +41 + (678)? — 678+1 
fi = 2924157727 


B=sx2+(s—-1I)x+sp-—1 


b = 678x(—80)2 + 677 x (-80) + 678*41 — 1 = 4312837 
Calculemos los dos valores para 
m=(sx2+s(s-2x+ps?-s+1)+T+sx2+(s-Dx+ps-—1 
my = 2924157727 + 12 + 4312837 = 35094205561 


n =(sêx2+s(s—-2)x+ps?-s+1)xT-—(sx2+(s— 1)x+ps) 
n, = 2924157727 + 12 — 4312838 = 35085579886 
Luego: 
Primer caso: 
bB=hATº+QB+HIT+AxX +x+p 


fp = 2924157727 + 122 + (2+ 4312837 + 1)+*12+(-80)2—- 80+41 
= 421182 227149 


Asi: 
(MN) +m+p=h+*b 
(35094205561)2 + 35094205561 + 41 = 2924157727 + 421182227149 
Su factorización completa es: 
(35094205561)? + 35094205561 + 41 = 563 * 748 103423 + 2924 157727 


Segundo caso: 


B=hATº-QB+IT+Ax +x+p 


fp = 2924157727 * 12? — (2+ 4312837 + 1)+12+(-80)2?—- 80 +41 
= 420975 210949 
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Así; 


(n)2+n,+p=h*f 
(35085579886)? + 35085579886 + 41 = 2924157727 * 420975210949 


Y esta es su factorización completa, o sea es un biprimo. 


8.3 Tercer Teorema de la factorización de Cordero en Z 


Este teorema se utiliza en la factorización de los números de la forma 4nº + 4p — 1 
con p € (2,3,5,11,17,41). Al igual que el primer y segundo teoremas de la factorización de 
Cordero en el conjunto de los números enteros, su factorización es de dos factores, y estos 
pueden ser compuestos, uno primo y el otro compuesto o los dos factores pueden dar 
números primos. 

Sea sx,TEZ,s*0, p E (2,3,5,11,17,41) 

1) Sin= foxtpot entonces 4n? + 4p — 1 es un número compuesto y dos de sus 
factores tienen la forma: 


f=x-3x+p+2 y b=x2+x+p 


2)Sin=[stx2+s(s—-2x+ps?-s+1]+T+ 


s(3s-2)x2+(3s2-8s+2)x+ps(35—2)-35+4 , 
Do R  |entonces4n?+4p—1 esun número compuesto 


y dos de sus factores tienen la forma: 


fi=sêx+s(s—-2)x+ps?-s+1 y 


fa 
=4ATº +8TB+(3s-— 2)2x? + (35 — 2)(3s — 8B)x 
' (3s2?—-8s+2)2+(3s— 2)[3s*(p — 1) — s2(2p — 11) — 135 + 2] 
s2 
mund p = [Hom ti tê mt sore tia] 


8.3.1  Aplicaciones del Tercer Teorema de la Factorización de Cordero 
en el Conjunto de los números enteros. 


Aplicación 1 
Sea x=90.p = 17, n = Erttet SÓSIA 4013 
f=x2-3x+p+2=902-3+90+17+2=7849 
b=x2+x+p=902+90+ 17 = 8207 


Luego: 4n? + 67 = 4+40132 + 67 = 7849 + 8207 
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Aplicación 2 
Seas=54x=-24p=11,T=13 


fh=sêx+s(s-2x+ps?-s+1 
=2916+*576+54+x52+x(-24)+11+2916-54+41 = 1644247 


 S0s=2)2º + (38º = 8842) + ps(3s=2)=35 +4 
E E BN. 2 A 
— 54+ 160+576 + 8318 + (—24) + 114 54+ 160 — 162 +4 


2 
= 2435945 


fo 
=4ATº+8TB+(3s— 2)2x? + (3s — 2)(3s — 8)x 
4 (3s?—-8s+2)2+(3s— 2)[3s*(p — 1) — s2(2p — 11) — 135 + 2] 
2 
=4+ 1644247 * 169 + 8x 13 + 2435945 + 1602 * 576 + 160 * 154 + (—24) 
a 69189124 + 160 + (472392 +10 — 2916 +11 — 13 +54+2) 


2916 


fp = 4+ 1644247 * 169 + 8+ 13 + 2435945 + 1602 * 576 + 160 + 154+(—24) 
+ 281129 


fp = 1379284621 
ó 


f = 4+ 1644247 * 169 — 8+ 13 + 2435945 + 1602 + 576 + 160 * 154 + (—24) 
+ 2811 = 872608061 


Luego: 

n=[s2x2+s(s-2x+ps?-s+1]+T 
N s(3s— 2)x2 +(3s2 —- 8s+2)x +ps(3s— 2) —-3s +4 
o Z 


n = 1644247 + 13 + 2435945 = 23811156 


n = 1644247 + 13 + 2435945 = 23811156 


Asi: 
4n?+4p—1 
4+ 238111562 + 43 = 1644247 + 1379284621 
[o 
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4 + 189392662 + 43 = 1644247 x 872608061 


8.4 Cuarto Teorema de la factorización de Cordero en Z 


Factorizar números enteros en sus factores primos utilizando fórmulas matemáticas 
no ha sido posible hasta el momento. Las fórmulas aquí publicadas permiten factorizar en 
dos factores algunos números enteros que tienen forma polinomial. El Cuarto Teorema de la 
Factorización de Cordero al igual que los anteriores permite factorizar en dos factores que 
pueden ser: primos los dos factores, uno primo y el otro compuesto o los dos compuestos. 

Seas xTEZSsHO ypeE(3,5,11,29) 

nsi n=(4s82x? +45x+2ps?+1)+T+ (25xº +x+sp) 


entonces: 2nº + p es compuesto y dos de sus factores tienen la forma: 
fj=4six2+4sx+2ps?+1 
2 
b= = = 26T2 +4BT+2x2+p 
1 


Con B=2sxº+x+sp 
2)Sin=(2x2+p)+*T+ (2sx2+x+sp) 
entonces: 2nº + p es compuesto y dos de sus factores tienen la forma: 


f=2x+p 
2 
b= TP =2fT2+4pT+2hs?+45x+1 
1 
Con B=2sxº+x+sp 


8.4.1 Aplicaciones del Cuarto Teorema de la Factorización de Cordero 
en los números enteros. 


Aplicación 1. 
Sea s=48x=23,p=29,T=12 


fj=4s2x?—-4sx+2ps?+1 


fi=4+(48)2(23)? —-4+48+23+2+29+482+1 = 5004481 
n= (48x? -4sx+2ps2+1)*T+ (25x? —-x+sp) 


n = 5004481 + 12 + (2+48+232 — 23 + 48+29) = 60105925 


n = 5004481 + 12 — (2+x48+232 — 23 + 48+29) = 60001619 
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2 
fo= = 2pT2 LAT + 2 +p 
q: 


f = 2+5004481+144+4+52153+12+2+232 +29 = 1443794959 


f = 2+*5004481+144-—4+52153+124+2+232 +29 = 1438788271 


Luego: 


2 * 601059252 + 29 = 5004481 + 1443794959 


2 * 600016192 + 29 = 5004481 * 1438788271 
Aplicación 2. 
Seas=8x=5,p=29,T7=11 
fi=2xº+p=2+25+29=79 
bB=2sx+x+sp=2+484*25+5+8+29=637 
n=(2x2+p)*T+ (2sx2+x+sp)=79x11+637= 1506 


n=(2x2+p)*T— (2sx2+x+sp)=79*11-—637=232 
Luego: 
b=2ATº+4BT+2fs2+4sx+1= 


2*79+x121+4+637+114+2+x79+64+4+8+5+1 =57419 


b=2fATº-4pT+2fs2+45x+1= 
2*79x121-4+637x11+2+79+64+4+x8+x5+1=1363 


Tenemos: 


2 + 1506” + 29 = 79+ 57419 


2+2322+29=79+136 


8.5 Quinto Teorema de la Factorización de Cordero en Z 


Sea sx,T,;reZ, sr+£0 ypeE(2,3,5,11,17,41) 
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Sin=(stx?+sx+p)(r*T+1)-r(sx+1)+1 entonces 
n2+(r—-2)n+pr? -r +1 es un número compuesto y dos de sus factores tienen 


la forma: 
fi=six2+sx+p 
b=A(r+aT+I2-r(2sx+ Dr4T+D+r? 
8.5.1  Aplicaciones del Quinto Teorema de la Factorización de Cordero 
en e! Conjunto de los números enteros. 


Aplicación 1 
Seas=-8, r=11,x=6 T=-4,p=41 


p(n)=n2+9n+4951 


n=(64+36-8+6+41)+*(-43)-11+(-8*6+D+1 


n=2297*(-43)4+11*47+1 
Donde: 


fi=2297 y n=-98253 
f = 2297x(—43)2 — 11443 +95 + 121 = 4202339 
Entonces: 
(-98253)? + 9 x (-98253) + 4951 = 9652772683 = 2297 * 4202339 


Y ambos factores son números primos. 
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RESUMEN-—Este informe surge por 
la necesidad de generar mejoras en la 
empresa, es un cambio necesario tanto en 
los procesos de realización del servicio, 
así como en la manera de capacitar a los 
empleados. Basado en la experiencia 
y el tiempo de trabajo en el ámbito 
correspondiente al giro de esta, tomaremos 
en cuenta los factores que influyen en 
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su desarrollo, dadas las premisas de la 
organización, que tiene una dirección clara 
a la productividad y la efectividad en los 
procesos. 

PALABRAS CLAVE-—Procesos, 
requerimientos, productividad, capacitación 
y mejora continua. 


INTRODUCCIÓN 


En México existe una gran 
diversidad de empresas, muchas de ellas 
carecen de algo tan simple como lo es la 
estructura organizacional, que es la parte 
medular para el crecimiento y expansión 
de esta. 

En mi experiencia en el campo 
laboral me he encontrado con muchas 
empresas que tienen este problema, y 
algunas se sostienen en el mercado, así 
como otras desaparecen, ya que carecen 
de un objetivo claro por el cual se guíen 
con convicción de cumplir las metas, no 
tienen reglamentos y/o procedimientos. 

Algunos empresarios no le toman 
importancia a este tipo de documentación 
que se requiere para elaborar una base 
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sólida como lo es la misión, visión, objetivos, reglamentos, políticas, los manuales de 
proceso y de bienvenida. 

El tiempo invertido en los procesos que se realizan en la empresa es más de lo 
necesario, dado que influyen varios factores en la actividad como la falta de sincronización 
entre los empleados, la falta de organización con los documentos que conforman los 
expedientes de los clientes, formatos incompletos con los cuales se solicita la información 
adecuada al cliente para aplicar a los programas SENTRI o GLOBAL ENTRY dicha 
información tiene que estar completa para introducirla en el sistema GOES (Global Online 
Enrollment Sistema), pero con el procedimiento establecido actualmente en la empresa hay 
mucho tiempo muerto y desperdicio de recursos. 

Hacer que la empresa funcione adecuadamente a través de los procesos que 
quedaran establecidos en el manual de optimización, aunado a esto el proceso de archivar 
expedientes, entre otros documentos, sea el que previamente se establece en el mismo 
manual. Así como también establecer el concepto de la empresa en el manual de bienvenida 
ante todos los empleados, son los objetivos de esta investigación. 

Al ver la situación en la que se encuentra la empresa es indispensable tomar acción 
en cuanto a la estructura organizacional y también establecer estándares por medio de los 
cuales la empresa se rija y se guíen los empleados para que no exista problema con la 
captura de información de los clientes. 

“Por lo que una organización puede estructurarse de diferentes maneras y estilos, 
dependiendo de sus objetivos, el entorno y los medios disponibles. La estructura de una 
organización determinará los modos en los que opera en el mercado y los objetivos que 
podrá alcanzar.” 


DESCRIPCIÓN DEL MÉTODO 


Los métodos para aplicar son: 

Descriptivo, porque mediante este tipo de investigación se utiliza el método de 
análisis y logra caracterizar el objeto de estudio. 

Cuantitativamente para datos numéricos, investigación empírica sistemática de los 
fenómenos sociales a través de técnicas estadísticas, matemáticas o informáticas. 

Cualitativamente porque es un método de investigación empleado en algunas 
disciplinas académicas, tradicionalmente en las ciencias sociales, sino también en la 
investigación de mercados y contextos posteriores. 

La recolección de datos se obtuvo mediante las etapas del proceso para la 
elaboración del manual, a través de la toma de tiempos al empleado, especificamente el 
tiempo que le tomaba realizar las solicitudes, dependiendo el grado de complejidad de esta, 
ya que en algunas fue necesario redactar un perdón, o la aclaración de alguna situación 
que afeció en la elegibilidad del aplicante. 


Estudios teórico-metodológicos en ciencias exactas, tecnológicas y de la tierra 2 Capítulo 3 


34 


Por otra parte, también se llevó a cabo el monitoreo de la competencia, en base a 
las promociones que publique, y sobre todo el tiempo que se tarda en realizar el trâmite 
completo. 

Por último, se midió la satisfacción del cliente, aplicando una pequefia encuesta, 
para que evalúe el servicio que se le brindo. 


MARCO TEÓRICO 


“La definición de la misión ayuda a clarificar las posiciones y creencias de los 
integrantes de la dirección de la empresa y a unificar criterios básicos. No es estrictamente 
necesario definirla para arribar a las políticas, especialmente en el caso de un único director 
o duefio, pero siempre es recomendable hacerlo, sobre todo para transmitir al personal los 
valores que constituyen la base de la cultura de la empresa. La visión es imprescindible, 
sobre ella se construyen las políticas de la empresa. Ejemplos de visión son: si la empresa 
es mediana, é seguirá siendo mediana o la dirección quiere construir una gran empresa”, si 
la empresa es líder en su país, é se plantea extender su liderazgo al continente o al mundo? 
Quizás la empresa nunca podrá realizar su visión, pero trabajará siempre en pos de ella. Si 
la política es el faro, la visión es la estrella.” (Nieves Cruz Felipe 2006) 

“Los valores, son principios éticos sobre los que se asienta la cultura de nuestra 
empresa y nos permiten crear nuestras pautas de comportamiento. 

No olvidemos que los valores son la personalidad de nuestra empresa y no pueden 
convertirse en una expresión de deseos de los dirigentes, sino que tienen que plasmar la 
realidad.” (Roberto Espinoza, 2015)2 

Un manual de procedimientos es un instrumento administrativo que apoya el 
quehacer cotidiano de las diferentes áreas de una empresa. 

En los manuales de procedimientos son consignados, metódicamente tanto las 
acciones como las operaciones que deben seguirse para llevar a cabo las funciones 
generales de la empresa. Además, con los manuales puede hacerse un seguimiento 
adecuado y secuencial de las actividades anteriormente programadas en orden lógico y 
en un tiempo definido. Los procedimientos, en cambio, son una sucesión cronológica y 
secuencial de un conjunto de labores concatenadas que constituyen la manera de efectuar 
un trabajo dentro de un âmbito predeterminado de aplicación. 

Todo procedimiento implica, además de las actividades y las tareas del personal, la 
determinación de los tiempos de realización, el uso de recursos materiales, tecnológicos y 
financieros, la aplicación de métodos de trabajo y de control para lograr un eficiente y eficaz 
desarrollo en las diferentes operaciones de una empresa. 

Las ventajas de contar con manuales de procedimientos son: 


*— Auxilian en el adiestramiento y capacitación del personal. 
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*— Auxilian en la inducción al puesto. 

* | Describen en forma detallada las actividades de cada puesto. 
* Facilitan la interacción de las distintas áreas de la empresa. 

*— Indican las interrelaciones con otras áreas de trabajo. 


* Permiten que el personal operativo conozca los diversos pasos que se siguen 
para el desarrollo de las actividades de rutina. 


* | Permiten una adecuada coordinación de actividades a través de un flujo eficien- 
te de la información. 


* Proporcionan la descripción de cada una de sus funciones al personal. 
* | Proporcionan una visión integral de la empresa al personal. 


* | Se establecen como referencia documental para precisar las fallas, omisiones y 
desempefios de los empleados involucrados en un determinado procedimiento. 


* | Son guías del trabajo a ejecutar3 


Las pequefias y medianas empresas, Pymes, tienen particular importancia para 
las economías nacionales, no solo por sus aportaciones a la producción y distribución de 
bienes y servicios, sino también por la flexibilidad de adaptarse a los cambios tecnológicos 
y gran potencial de generación de empleos. Representan un excelente medio para impulsar 
el desarrollo económico y una mejor distribución de la riqueza. Hoy día, los gobiernos 
de países en desarrollo reconocen la importancia de las Pymes por su contribución al 
crecimiento económico, a la generación de empleo, así como al desarrollo regional y local. 
Las Pymes en México constituyen el 97% del total de las empresas, generadoras de empleo 
del 79% de la población y generan ingresos equivalentes al 23% del Producto Interno Bruto 
(PIB), lo anterior es una clara serial de debemos poner atención a este tipo de empresas y 
verlas como lo que en realidad son: la base de la economía mexicana. (INEGI, 2005). 


MARCO CONTEXTUAL 


Global Express Access es una empresa dedicada al servicio en trâmites migratorios, 
como lo son la visa, residencias y se especializa más en la línea express, tiene laborando 
en el mercado nueve afÃios que respaldan a la empresa con la experiencia adecuada y 
necesaria para llevar a cabo un excelente servicio y atención al cliente. 

Los programas Sentri y Global Entry, que se encuentran en el GOES (Global Online 
Enrollment System) son programas para viajeros confiables a lo largo del mundo, además 
de los Estados Unidos y Canadá, se han integrado al programa otros países como Holanda 
y Corea. Este programa permite que además de cruzar por los puentes internacionales, 
al internarte en los aeropuertos de alguno de estos países también se pueda evitar la 
inspección, evitando así retrasos y perdidas de vuelos de conexión. Muchos de estos 
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aeropuertos han sido previamente seleccionados por el volumen de personas que viajan 
fuera de Estados Unidos, en dichos aeropuertos encontraras una serie de kioscos que te 
permiten emitir un ticket para que pueda ser fácilmente transferido al área de recepción de 
maletas. Además, en algunas ocasiones hasta evitar las inspecciones antiterroristas que 
te pueden resultar en que pierdas tus vuelos de conexión. Estados Unidos está ofreciendo 
una serie de programas de viajeros confiables para ciudadanos mexicanos y canadienses. 

Global Express Access ofrece asistencia personalizada y se le da atención especial 
a cada uno de los clientes, asesorando paso a paso al usuario, desde que comienza la 
inscripción o renovación hasta que la tarjeta les es entregada, de igual manera cualquier 
duda, aclaración o cambio se efectúa en tiempo y forma. 


MARCO CONCEPTUAL 


SENTRI: conocido también como línea express permite a los viajeros que crucen 
la frontera de Estados Unidos vía terrestre desde México o Canadá por cualquiera de los 
puertos de entrada que cuenten con una línea dedicada para este servicio cruzar de una 
manera rápida y expedita con el uso de su tarjeta SENTRI. Este programa debe de ser 
utilizado inicialmente por los usuarios que no hayan cometido ninguna falta o delito mayor 
en ninguno de los países no solo en Estados Unidos 

GLOBAL ENTRY: es un programa para viajeros confiables a lo largo del mundo. Este 
programa permite que además de cruzar por los puentes internacionales, al internarte en 
los aeropuertos de alguno de los países registrados también puedas evitar la inspección, 
evitando así retrasos y perdidas de vuelos de conexión. 

CBP: la Oficina de Aduanas y Protección Fronteriza de los Estados Unidos (U.S. 
Customs and Border Protection, CBP) es una agencia del Gobierno de los Estados Unidos. 
Es la que controla la protección de los puertos internacionales de cruce en USA realiza una 
estricta revisión de registros personales y califica a los usuarios como viajeros confiables. 
La CBP protege contra narcotráfico, terrorismo, y tráfico ilegal de personas. También tiene 
competencia en el comercio e inmigración legítimos. Tiene más de 52.000 empleados, 
incluyendo más de 22.000 funcionarios y especialistas en agricultura, más de 17.000 tropas 
de frontera, y 1.000 agentes del aire y agentes marinos. Tiene su sede en Washington, DC. 


RESULTADOS, CONCLUCIONES Y RECOMENDACIONES 


Conclusiones 


1: Los resultados arrojados muestran que se cuenta con buen servicio 
considerando el tiempo de proceso con el trámite de la línea express, es 
regular considerando todos los factores de espera, pero lo ideal sería causar 
una buena impresión con un excelente servicio. 
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2: Muchos de los factores que influyeron en el tiempo de proceso fue la falta 
de organización de la empresa ya que se perdía rápidamente el seguimiento 
de los casos. 


3: La capacitación a los empleados era deficiente ya que se le dedicaba muy 
poco tiempo a ello y el empleado aprendia con la práctica, mismo detalle que 
conileva mucho tiempo. 


4: En base al análisis de los resultados los clientes consideraron que el tiempo 
que se invierte en la gestión de inscripción y/o renovación fue en su mayor 
parte “tardado” y afecta los intereses y la imagen de la empresa, ya que la 
publicidad que hace un cliente satisfecho es la mejor. 


5: La inducción de los nuevos procesos se reflejó considerablemente en el 
flujo de aplicaciones certificadas a la semana. 


6: La definición de los puestos y de tareas favoreció la convivencia entre los 
empleados ya que no hubo confusiones ni trabajos duplicados. 


7: Como se cuenta con un proceso ya establecido de como elaborar el servicio 
es benéfico ya que agiliza su elaboración, lo que reduce posibles problemas 
que puedan llegar a surgir. 


8: Muchos de los clientes les son leales a la empresa, lo que es bueno ya que 
muchos de ellos recomiendan el servicio, sin embargo, ellos reciben un trato 
especial por ya permanecer afios en la lista de clientes activos. 


9: A pesar de las circunstancias hubo mejoras en la empresa en general, que 
con el paso del tiempo se podrán reforzar. 


10: Un punto muy importante en cada empresa son los empleados, ya que 
son una vértebra muy importante, es indispensable que los clientes se vayan 
satisfechos con el servicio, pero también es necesario que los empleados se 
encuentren motivados en su trabajo. 


Recomendaciones 


1: Se recomienda un buzón interno para atender las peticiones o los puntos de 
vista de los empleados ya que ellos son los que están es servicio directo con 
el cliente y de ante mano saben las necesidades de las personas que atienden. 


2: Es importante que los puntos de vista expuestos los analicen y sean 
tomados en cuenta y verifiquen si se pueden hacer cambios para mejorar el 
servicio. 


3: Tomando en cuenta que, si les afecta hasta cierto punto los nuevos 
procesos al elaborar las aplicaciones, podrán tomar un tiempo designado 
exclusivamente para realizar prácticas de prueba y error. 


4: La nueva designación de tareas y actividades les ahorrará tiempo y se podrá 
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invertir en nuevos proyectos para la empresa. 


5: En ocasiones será necesario realizar un servicio express ya que es por 
petición del cliente se recomienda estar muy activo para brindar el mejor 
servicio. 


6: Seguir el reglamento establecido para poder impulsar el orden en cuanto a 
los empleados y el ambiente sea más próspero. 


7: Es de suma importancia que estén bien establecidos los lineamientos 
sobre el servicio, se les recomienda seguir los estándares de los procesos 
ya establecidos y apoyar a los empleados con los mismos para no olvidar los 
puntos importantes, al momento de Ilevarlos a cabo con el cliente. 


8: Para evitar contradicciones con los empleados es importante no perder de 
vista los objetivos y el reparto de las actividades correspondientes a cada uno. 


9: Tomando en cuenta que se asesora al empleado para que siga los 
lineamientos al pie de la letra para que no le afecte las amonestaciones, ni en 
su sueldo ni en su récord como empleado. 


10: Partiendo de las amonestaciones, también premiar y reconocer la labor de 
los empleados en la empresa, ya que, sin duda, son el alma de la misma. 
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CAPÍTULO 4 


METODOLOGIA SIMPLIFICADA PARA 
AVALIAÇÃO GEOTÉCNICA DE UM DEPÓSITO 
CONTROLADO DE ESTÉRIL (DCE) 


Fernando Alves Cantini Cardozo 


RESUMO: (0) trabalho apresenta 
uma | metodologia simplificada para 
desenvolviemnto de estudo geotécnico e 
avaliação de estabilidade de um Depósito 
Controlado de Estéril (DCE) de uma Mina 
de calcário na cidade de Dom Feliciano/ 
RS. Com base em revisão bibliográfica 
são feitas considerações a respeito de 
propriedades geotécnicas de materiais e 
de geometrias adotas. Tais parametros 
embasaram a avaliação de estabilidade 
feita sobre perfis realistas da DCE e sobre 
uma proposta de ampliação. As conclusões 
quanto a estabilidade são confrontadas e 
corroboradas por metodologia quantitativa 
de classificação de DCE. Os resultados 
obtidos corroboram pela viabilidade da 
adoção de uma geometria otimizada, com 
altura de 25 metros para o primeiro patamar 
e 15 metros para o segundo. Conclui-se 
que a metodologia e sequencia de passos 
adotados na avaliação pode ser replicada 
para outras DCE e até mesmo servir de 
indicativo para estruturas semelhantes. 

PALAVRAS-CHAVE: Metodologia, Estéril 
de Mina, Avaliação Geotécnica, Análise de 
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Estabilidade 


INTRODUÇÃO 


Pilhas de Estéril figuram juntamente 
com Barragens de Rejeito entre as maiores 
estruturas artificilamente contruidas. 
Tais estruturas podem facilmente atingir 
centenas de metros de altura a dependenr 
da escala de produção adoptada na 
mineração. Naturalmente é consituida 
de rejeito proveniente da lavra, sendo 
principalmente da decapagem e por vezes 
de processos de beneficialemento. Diferem- 
se principamente dos rejeitos destinados 
a barragens pelo baixo teor de água e 
viabilidade pelo empilhamento. Assim 
como demais estruturas de uma mina, os 
Dépositos de Estéril, podem apresentar 
riscos geotecnicos e dependendo de suas 
caracterisitica podem vir a necessitar de 
extenssos estudos para sua adequada 
avaliação. Wang et al. (2020)e Ulusay et 
al. (1995) apresentam exemplo e outros 
casos onde a rupturas de DCE mobilizaram 
milhares de metros cubicos de materil. 


Neste estudo é apresentada a 
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proposta de avaliação geotécnica feita para o depósito controlado de estéril (DCE) da Mina 
Gaspar Simões, localizada no município de Dom Feliciano/RS, teve como objetivo avaliar a 
estabilidade geotécnica da estrutura citada na situação atual e considerando sua proposta 
de ampliação de capacidade. O DCE fica localizado na área útil do empreendimento e é 
constituido do estéril da lavra, constituido predomiantemente de rochas encaixantes ao 
minerio marmore calcitico; sendo, granitos, gnaisses, mamores e saprolito. Nos estudos 
anteriormente apresentados para o licenciamento ambiental do empreendimento, foi 
proposta de forma empírica, uma geometria conservadora, e que considerou a época uma 
reconfiguração do DCE. A proposta original, previa uma configuração de com taludes de 10 
metros de alturas, inclinação de 45º e bermas de 4 metros. Esta proposta inicial demandaria 
uma extensa intervenção na configuração atual e extensos custos decorrestes deste. 

Assim, este estudo visa certificar a estabilidade da DCE nas condições atuais de 
geometria e verificar a viabilidade de execução de sua ampliação e configuração final, 
considerando a necessidade de avanço de lavra. 
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Figura 1 — Locaçização da Pilha de Estril da Mina Gaspar Simões 


METODOLOGIA 


Para o desenvolvimento da Avaliação da situação atual do DCE e proposição de 
configuração final, optou-se por seguir a metologia desenvolvida exposta no fluxograma 
da Figura 1. A metodologia usa em paralelo a avaliação quatitativa, com determinação 
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de Fator de Segurança; e Avaliação qualitativa, está com intuito de classificar o DCE e 
subsidiar o projeto. Para a valiação quantitativa optou-se pela analise por equilibrio limite 
para calculo do Fator de Segurança. Já a análise qualitativa consistiu na aplicação da 
metodologia adaptada por Aragão (2008). 

A metodoliaga em parte se baseia na definição de perfis representativos e 
avaliação de modelos geotecnicos baseados nestes. Destacando que a utilização de 
Modelos digitais de Terreno permite extração de perfis com grande detalhamento, e vem 
sendo empregada em diversos trabalhos de geotecnia. Na metodologia em questão, o 
levantamento com VANT e geração de MDE também subsidia avaliação qualitativa, ao 
fornecer um maior detalhamento e confiabilidade a interpretação das feições da estrutura. 
Assim, complementando as inspeções de campo, implicitas as etapas de identificação e de 
mapeamento da área. Como refencia ao uso de VANT e metodologia de mapeamento com 
aplicação em geotecnia podemos citar os trabalhos de Zekkos et al. (2018), Fu-Hsuan et 
al. (2018). 

Por fim, os resultados da avaliação devem ser analisados e definidos, se necessario, 
aprofundamentos de estudos e/ou alterações de projetos, como alteração de geometria. 
Destaca-se que a analise qualitativa, conjuntamente com o FS obtido, permite validar 
sobretudo os parametros utilizados nos modelos. Uma vez, que pode se entender pela 
necessidade de aprofundamentos de estudos e ou instrumentações. Na Fluxograma da 
metodologia aplicada, as setas tracejadas indicam fluxos opcionais. 
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Identificação do DCE 


Avaliação da documentação 
disponível 
*Estudos anteriores, cartografia, bibliografia 


Levantamento Aéreo - Avaliação Qualitativa 
Modelagem Digital do Terreno Sistema de Classificação Aragão (2008) 


para avaliação geotécnica di 
Definição da Classe do DCE 
Risco e ruptura/Recomendações 
Revisão bibliografica e criação dos 
modelos geotécnicos 


Avaliação dos resultados e reavaliação do projeto 
(aprofundamento de estudos e/ou redefinição do projeto) 


Figura 2 - Metodologia proposta para avaliação do DCE 


Situação da Pilha e Estudos Anteriores 


A Pilha de Estéril da Mina Gaspar Simões encontra-se em operação, tendo sido 
reliazados estudos para fins de licenciamento ambiental da atividade no ano de 2021. 
Tendo sido realizados sondagem e proposição de reconficguração para uma geometria 
mais conservadora. Destaca-se que embora tenha sido realizado ensaio de campo (SPT) 
e utilizado para inferência de paramêtros, este atingiu a profundidade de 0,80 metros na 
DCE; atingindo como esperado o impenetravel — uma vez que o depósito controlado de 
estéril (DCE) é composto de blocos de tamanhos variados e de solo compactado. Como 
proposta de configuração original, o projeto apresentava taludes de 45º e altura máxima 
de 10 metros subdivididos por bermas de 4 metros. Assim, seria necessaria uma extensa 
reformulação nos perfis do DCE atual com alturas entre 19 e 25 metros. Em tal proposta 
o FS chegaria a 3.29. Já em 2022 foi realizado nova avaliação com finalidade verificar a 
estabilidade da situação atual e de nova proprosição de perfil final. Como base do estudo, 
fora feita das documentações disponives, caminhamento pela estrutura e mapeamento com 
VANT. Posteriormente as imagens foram processadas em ambiente GIS para verificação 
dos perfis e avaliação das geometrias. A Figura 2 apresenta o MDE da Pilha e os principais 
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perfis identificados. 

A DCE apresenta faces com inclinação limitada a 35º, sendo no setor leste de tem 
altura máxima de 25 e média de 20 metros, com inclinação na ordem de 30º. No setor oeste, 
a DCE apresenta ângulo de face na ordem de 34º a 38º, sendo sua altura de bancada 
entre 12-16 metros. Considerando a face da DCE, ela apresenta patamar em formação 
ainda significativamente afastado da borda, não influindo na avaliação da situação atual. 
Destaca-se que para o setor oeste é projetada ampliação, a qual em sua configuração final 
pretende regularizar a inclinação máxima em até 35º. 
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Figura 3 - Vista em planta contendo a localização dos perfis avaliados da DCE. 


Com base no Modelo Digital de Elevação (MDE), pode-se estimar a declividade 
primitiva da área da DCE. Na área de interesse verifica-se uma declividade leve a moderada, 
na ordem de 10% na porção leste da DCE (área dos perfis 1,2 e 3) e da ordem de 5-7% na 
porção oeste (área dos perfis 4, 5 e 6). Não se observam surgências ou acúmulos de água 
ou umidade no entorno da DCE que permitam supor qualquer grau de saturação interna. 
Tais informação são relevantes para a avaliação qualitativa e quantitativa. 


AVALIAÇÃO GEOTECNICA 

O estéril de mineração por natureza é constituído da decapagem do minério sendo 
constituído majoritariamente por fragmentos rochosos em meio a matriz de saprólito. É 
característica típica a alta granulometria do estéril, devido a economicidade aplicada no 
desmonte do material da decapagem, chegando a granulometria de pedregulho, matacos e 
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blocos. Dada a alta granulometria é comum a consideração para projetos o próprio ângulo 
de repouso dos materiais, uma vez que a deposição se dá pelo lançamento do estéril. 
Todavia, para avaliação geotécnica de estabilidade é pertinente a avaliação dos demais 
parâmetros como densidade natural, coesão e ângulo de atrito. 

Diversas metodologias são apresentadas na literatura, para avaliação de estabilidade 
de depósitos de estéril. Normalmente esta pode ser realizada por métodos numéricos (de 
elementos finitos, discretos, entre outros) e principalmente por meio do método de equilíbrio 
limite, com o critério de falha de Mohr — Coulomb, por exemplo, considerando ruptura por 
cisalhamento interno do material, em uma superfície circular preferencial. 

Referente aos parâmetros do modelo, depósitos de estéril, apresentam grandes 
limitações quanto a determinação/estimativa direta de suas propriedades de resistência, 
sendo limitados o acesso a dados e publicações quanto o resultado de ensaios de campo e/ 
ou de laboratório. Portanto, são comuns estimativas de resistência baseadas no ângulo de 
repouso observado e ou dados teóricos por bibliografia, para materiais de granulometrias 
semelhantes. 

Cabe salientar que ocorre divergência quanto ao valor considerado para a coesão, 
cuja ocorrência se relaciona com características da matriz onde se encaixam os fragmentos 
maiores; e/ou a escala das estruturas. A Tabela 1 abaixo apresenta levantamento de 
parâmetros de resistência considerados em bibliografia pelo critério de Mohr-Coulomb. 


Autor especíto (av | Sogsão | Ângulo de 

Saliba et al. (2016) 22 8 28 
Gomez, Garcia e Romanel (2016) - 0 34 
Castro (2022) 19 10 32 
Zou et al. (2018) 22,1 35 35 
Stormont e Farfan (2005)* 16,5 4,8 37 
Huang (1983) apud Stormont e Farfan (2005) - - + de 34 
Das (2002) apud Stormont e Farfan (2005) - - 34-48 
McDaniel and Blair, (1987) apud Stormont e Farfan & - 34 
(2005) 
Mariani (2017) 26,29 18,7 37,4 

Cruz (1996)** - 0-2 33-55 


*Material estéril de rocha calcária. **Enrocamento em geral. 


Tabela 1 - Parametros de resistencia para estéril de mina enrocamento. 


Destaca-se os parâmetros médio de ângulo de atrito superior a 35º e sendo 
recorrente o valor de 37º. Quanto a coesão, está apresenta maior variação sendo nula 
nos materiais puramente granulares, mas usualmente referenciada com o valor de 10 kPa 


nos trabalhos que consideram variação de granulometria, com destaque para o trabalho 
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de Stormont e Farfan (2005); também de Mariani (2017) que avaliou material visualmente 
parecido com a composição da matriz do estéril. 

Quanto ao peso específico ou densidade natural, se observa uma tendência a 
escala da ordem de 20 kN/mº3, valor que é condizente com densidade de rochas (25-27 kN/ 
mº) combinadas com espaços vazios e solo/rocha intemperizada preenchendo os espaços 
vazios. Vale destacar que o estudo anterior considerado, com estimativa de parâmetros 
a partir do valor de SPT, utilizou o valor de 75 kPa de coesão e 40,7º de ângulo de atrito 
interno; onde se observa que principalmente a coesão, apresenta um valor extremamente 
alto. Sendo assim, a reconsideração dos parâmetros geotécnicos da DCE, para a avaliação 
de estabilidade foi feita de forma a aproximá-los da realidade, tornando a avaliação atual 
mais conservadora em relação à anterior, porém mais próxima da realidade. Foi considerado 
Peso específico de 20 kN/m?, Coesão de 10 kPa e Ângulo de atrito: 37º 


Análise de Estabilidade 


Para as análises, utilizou-se o software Slide 6, estruturando modelos geotécnicos 
que evitassem efeitos de borda. Para a análise de estabilidade, definiu-se a subdivisão das 
superfícies de ruptura em 25 lamelas. Para modelagem da linha piezométrica, considerou- 
se duas hipóteses conforme preconização da NBR 13029, com cenário real (seco/drenado) 
sem ocorrência de nível d'água na DCE e como cenário hipotético com ocorrência de nível 
d'agua. Como critério de ruptura, utilizou-se o critério de Mohr-Coulomb. 

Os fatores de segurança foram calculados pelos Métodos de Bishop Simplificado, 
considerando superfícies circulares, considerando ser o modelo de ruptura de larga utilização 
na bibliografia consultada. Como critério de linha freática crítica, considerou-se saturação 
a meia altura de aterro, destacando-se ser uma avaliação conservadora, uma vez que, sua 
ocorrência é improvável dada a drenagem externa da DCE e da mesma ser parcialmente 
alto drenante dada a variação granulométrica de sua constituição. Assim afasta-se a 
hipótese de saturação na DCE, embora por normativa esta deva ser considerada em caso 
hipotético. Nos resultados dos modelos, estão expostas a superfície de ruptura crítica e a 
linha piezométrica considerada para situação crítica. A seguir estão expostos os resultados, 
com as geometrias da PDE atual e projeção; e modelos geotécnicos considerados. 

Referente a recomendações de Fator de Segurança a NBR 13029:2017, apresenta, 
conforme Quadro 1 abaixo, os Fatores de Segurança recomendados para pilhas de estéril 
de mineração. A altura é limitada a 200 metros, não especificando valores de inclinação e/ 
ou alturas de bancadas. Todavia, destaca-se que estes devem ser observados quanto a FS 
da estrutura, operacionalidade e otimização; no intuito minimizar danos ambientais. 
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Ruptura do Talude Geral Superfície FS mínimo 
Freática normal 1,5 
Freática crítica 1,3 

Ruptura do Talude Geral Face FS mínimo 
Predominantemente de solo 1,5 
Predominantemente de rocha 1,3 


Quadro 1 - Fatores de segurança da ruptura do talude geral e do talude entre bermas (modificado de 
ABNT 13029:2017) 


Cabe destacar que Peck (1969), comenta que estudos geotécnicos podem ser 
divididos quanto ao detalhamento da investigação em três níveis: 
1- Baseados em investigações geotécnicas limitadas e adotando abordagens de 
projeto conservadoras e altos fatores de segurança; 


2- Baseados em investigações geotécnicas limitadas e avaliando recomendações 
baseadas em práticas regionais; 


3- Executando investigações detalhadas. 

Compreende-se que o estudo em questão se apoia nos níveis 1 e 2, apresentando 
investigações limitadas e práticas regionais - dadas as investigações feitas - e sendo 
baseados em práticas regionais (quanto a geometria e profundidades de solo), tanto na 
inferência de valores, quanto na observância da estabilidade da estrutura em questão já 
existente e de outras em situação semelhante. 


Resultados Situação 


Conforme perfis numerados na figura 2, extraídos do MDE. 


Figura 4 — Perfil 1 - situação normal Figura 5 — Perfil 1 - situação com freática crítica 
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Figura 8 — Perfil 3 - situação normal Figura 9 — Perfil 3 - situação com freática crítica 


Figura 12 — Perfil 5 - situação normal Figura 13 — Perfil 5 - situação com freática crítica 
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Figura 14 — Perfil 6 - situação normal Figura 15 — Perfil 6 - situação com freática crítica 


A tabela abaixo apresenta os fatores de segurança obtidos nas condições de 


avaliação. 
Perfil | Perfil 1 Perfil 2 Perfil 3 Perfil 4 Perfil 5 Perfil 6 
FS. 
Obtido 1.612 1.628 1.588 1.684 1.625 1.602 
Cenário Crítico | 1.385 1.412 1.488 1.655 1.432 1.406 


Tabela 2 - Resumo dos Fatores de Segurança para os Perfis Atuais 


Resultados para Ampliação e Configuração Final da PDE 


Para a ampliação do setor leste do DCE atual, foi considerada execução de patamar 
adicional com altura de 15 metros e inclinação de 35º, com berma mínima de 5 metros 
(espaço para eventuais manutenções), dividindo o talude frente ao primeiro patamar.A 
avaliação de estabilidade apresentou FS de 1,5 para talude global na situação normal e 
1,3 para situação hipotética com ocorrência de nível freático. As Figuras abaixo se referem 


aos perfis 1 e 2. 


Figura 16 — Perfil 1 Ampliação - situação normal Figura 17 — Perfil 1 Ampliação - situação crítica 
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Figura 18 — Perfil 2 Ampliação - situação normal Figura 19 — Perfil 2 Ampliação - situação crítica 


Para as áreas do DCE atual onde fora demarcado avanço vertical e para execução 
do novo DCE no setor norte foi considerado a configuração de primeiro patamar com talude 
de 25 metros e 35º de inclinação, e segundo patamar com talude de 15 metros e inclinação 
de 35º; utilizando berma de 5 metros. 


AVALIAÇÃO QUALITATIVA 


Para a avaliação de estabilidade do DCE em caráter qualitativo, baseado em 
informações visuais, feições geomorfológicas, geometrias do DCE e terreno, drenagem, 
hidrologia, hidrografia local e regional. Aragão (2008) e Aragão e Oliveira Filho (2011), 
apresentam metodologia baseada sistema de classificação de pilhas desenvolvido para 
o governo canadense com base na estabilidade física de uma pilha de estéril (BC MINE 
WASTE ROCK PILE RESEARCH COMMITTEE, 1991). Tal metodologia é recomendada 
para a avaliação da estabilidade do DCE, por ser útil e de fácil aplicação permitindo a 
classificação e enquadramento do DCE. Adicionalmente o sistema também apresenta 
conforme a classificação sugestão de medidas de monitoramento. As tabelas completas 
podem ser encontradas em Aragão (2008). 

As Tabelas a seguir apresentam a aplicação da metodologia para o DCE da Mina 
Gaspar Simões, considerando a situação atual, ampliação projetada e execução da nova 
pilha projetadas para a área norte. 


Estudios teórico-metodológicos en ciencias exactas, tecnológicas y de la tierra 2 Capítulo 4 


51 


Critério Condição/Descrição Pontuação 
Mina GS 
Altura <50m 0 
Volume da Pilha < 1*1016 m?/banco 0 
Inclinação do Talude 26º-35º 50 
Inclinação Fundação 10º 0 
Confinamento - Bancos ou terraços naturais nos taludes; 50 
- Taludes com inclinação uniforme, limitados por 
topografia natural diversificada; 
-Empilhamento de estéreis em encostas, em vales 
abertos, ou transversais a vales. 
Fundação - Materiais de Fundação tão ou mais resistentes que os | 0 
da pilha; 
- Não sujeita a efeitos adversos da poropressão; 
- Sem estruturas geológicas desfavoráveis. 
Qualidade do material da | - Moderadamente resistente, durabilidade variável; 100 
pilha - 10 a 25% de finos 
Método construtivo - Bancos ou camadas não muito espessos (<25m de 0 
espessura), plataformas largas; - Disposição ao longo 
das curvas de nível - Construção ascendente - Wrap- 
arounds ou terraços 
Taxa de disposição - < 25m3 /banco por metro linear de crista por dia; 0 
- Taxa de avanço da crista <0.1 m por dia 
Sismicidade Baixa 
Condições piezométricas | - Baixas pressões piezométricas, nenhuma surgência 0 
e climáticas na fundação; 
- Improvável desenvolvimento de superfície freática no 
interior da pilha; 
- Precipitação limitada; 
- Infiltração mínima dentro da pilha; 
- Nenhuma camada de neve ou gelo na pilha ou 
fundação 
Somatório 200 


Tabela 3 - Pontuação para DCE Mina Gaspar Simões 


Classificação 
da Ruptura 
estabilidade 


Potencial de 


Faixa de 
pontuação 


Nível de Esforço Recomendado 
para Investigação, Projeto e 
Construção 


Classificação da 
PDE 


| Desprezível 


- Reconhecimento básico do local, < 300 
documentação de referência; 

- Programa mínimo de ensaios de 
laboratório; 

- Rotina de checagem de 
estabilidade, possivelmente usando 
ábacos; 

- Restrições mínimas na 
construção; 

- Monitoramento apenas visual. 


Mina Gaspar 
Simões 


200 


Tabela 4 - Classificação da PDE 
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Classe da Monitoramento Descrição Frequência Relatório 
Pilha Instrumental 
| - Inspeção visual; Encarregado de 4 horas Relatório Diário 
- Piezômetros onde é | turno faz a inspeção 
aplicável. — — 
Inspeção periódica Anual Relatório Anual 
detalhada feita pelo 
Engenheiro 


Tabela 5 - Recomendações monitoramento PDE 


Verifica-se que o DCE atual se enquadra como Classe le risco de ruptura desprezível 
segundo a metodologia. Mesmo considerando as ampliações e configuração final do 
DCE, o mesmo mantem-se na Classe |, destacando que até Classe Il, o risco atribuído a 
potencial de ruptura é baixo. Sendo recomendável nestas Classes o monitoramento visual, 
avaliações quantificas por ábacos e instrumentações quando cabível. Mesmo na situação 
de ampliação e configuração final o DCE tende a manter seu enquadramento, considerando 
a geometria. 


Configuração geométrica 

A tabela abaixo, apresenta a configuração geométrica proposta com base na 
avaliação realizada. Tal configuração considerou as verificações qualitativas e quantitativas, 
as quais indicam viabilidade pela execução nas configurações propostas. Destaca-se que 
a proposição da geometria além da viabilidade geotécnica, se sustenta na otimização de 
espaço e minimização de impactos, reduzindo movimentações de terra e impactos em 
área planificada; uma vez que se entende que na configuração proposta atualizada a 
estabilidade é garantida (se enquadrando em norma). 


DCE atual + ampliação Altura bancada Inclinação banca 
Primeira bancada Até 25 metros Até 35º 
Segunda bancada Até 15 metros Até 35º 


Tabela 6 - Configuração dos DCE. 


CONSIDERAÇÕES FINAIS 


A metodologia adotada se mostrou adequada para o caso em questão. Os resultados 
das avaliações realizadas corroboraram pelo entendimento de estabilidade do DCE na 
situação atual e de sua ampliação proposta. A avaliação de equilíbrio limite, feita sobre os 
perfis, apresentou FS condizentes com o sugerido na NBR 13029:2017. Tal entendimento 
é reforçado pela avaliação qualitativa, que em função de características da Pilha/Depósito 
apontou para sua estabilidade — classificando de baixo risco. Sendo sugerida inspeções 
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visuais rotineiras e inspeções de maior detalhamento anualmente. 

Considerando que avaliação geotécnica aponta para estabilidade e baixo risco da 
situação atual e da ampliação proposta, podemos concluir que a adoção de geometria 
mais conservador do projeto original se configuraria em um custo financeiro e ambiental 
desnecessário. Sendo que a configuração atual já cumpre as necessidades de estabilidade 
e operacionais. 

Referente a metodologia adotada, identifica-se potencial para aplicações na 
estabilidade de outros depósitos ou mesmo outras situações, como em taludes escavados 
e barragens. Desde que observadas as necessárias adaptações quanto a avaliação 
qualitativa. Ao exemplo deste estudo, a avaliações geotécnica de estabilidade pode validar 
ou apontar para configurações que sejam mais econômicas ao empreendimento e menos 


impactantes ao meio ambiente, ao se evitar estruturas superdimensionadas. 
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RESUMEN: La zona de estudio corresponde 
a los relieves de lomadas constituidos por 
series alternantes de limolitas y areniscas 
finas, de la Formación El Paso (Pleistoceno), 
ubicados al sudoeste de la provincia de 
Salta, Argentina, entre la localidad de 
Cafayate y el paraje Lorohuasi. Los estudios 
geomorfológicos realizados permitieron 
determinar que la continuidad del primitivo 
depósito lacustre fue afectada por ciclos 
de humectación-desecación resultantes de 
las precipitaciones estacionales, propias 
del clima seco-árido que caracteriza a 
la zona desde el Holoceno al presente 
y a la intensa amplitud térmica diaria, 
entre 12ºC y 36ºC. Las precipitaciones 
estacionales favorecen el desarrollo de 
una intensa erosión e incisión pluvial, en 
los depósitos lacustres, generando típicos 
surcos de erosión conocidos como piping. 
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Los surcos de escurrimiento se originan 
cuando el agua erosiona y elimina el 
material que se encuentra en la parte 
superior del depósito. Estos sedimentos 
son transportados hacia los laterales y la 
parte baja formando pequefios abanicos. El 
material retransportado hacia los laterales 
forma sobre el depósito original una nueva 
cobertura de material pelítico. Las pelitas del 
depósito original se diferencian de las del 
depósito secundario, porque las primeras 
presentan laminación, mientras que las 
segundas tienen una estructura maciza. 
PALABRAS CLAVE: Sedimentos lacustres, 
Hidroclástismo, Surcos de erosión. 


RELIEF MODELING IN OUTCROPS 
OF THE EL PASO FORMATION, 
BETWEEN CAFAYATE AND 
LOROHUASI, SALTA, ARGENTINA 
ABSTRACT: The study area corresponds to 
the reliefs of hills assembled by alternating 
series of silistones and fine sandstones, 
of the El Paso Formation (Pleistocene), 
located to the southwest of the province 
of Salta, Argentina, between the town of 
Cafayate and the Lorohuasi area. The 
geomorphological study carried out made it 
possible to determine that the continuity of 
the primitive lacustrine deposit was affected 
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by cycles of humidification-desiccation resulting from seasonal rainfall, typical of the dry-arid 
climate that characterizes the area from the Holocene to the present and, to the intense 
daily thermal amplitude, between 12ºC and 36ºC. Seasonal rainfall favors the development of 
intense erosion and pluvial incision, in the lacustrine deposits, typical erosion grooves known 
as piping will result. Runoff grooves are created when water erodes and removes material 
on top of the reservoir. These sediments are transported to the small sides and the bottom 
forming fans. The material retransported to the sides forms a new cover of pelitic material over 
the original deposit. The pellets from the original deposit differ from those from the secondary 
deposit, because the former present lamination, while the latter have a solid structure. 
KEYWORDS: Lacustrine sediments, Hydroclasticism, Erosion furrows. 


MODELAGEM DE RELEVO EM AFLORAMENTOS DA FORMAÇÃO EL PASO, 
ENTRE CAFAYATE E LOROHUASI, SALTA, ARGENTINA 

RESUMO: A área de estudo corresponde aos relevos de colinas constituídos por séries 
alternadas de siltitos e arenitos finos, da Formação El Paso (Pleistoceno), localizados 
ao sudoeste da província de Salta, Argentina, entre a cidade de Cafayate e a área de 
Lorohuasi. Os estudos geomorfológicos realizados permitiram constatar que a continuidade 
do depósito lacustre primitivo foi afectada pelos ciclos de humidificação-dessecação 
resultantes das precipitações sazonais, típicas do clima árido-seco que caracteriza a zona 
desde o Holocénico até ao presente, e pela intensa amplitude térmica diariamente, entre 
12ºC e 36ºC. As chuvas sazonais favorecem o desenvolvimento de intensa erosão e incisão 
pluvial nos depósitos lacustres, gerando sulcos de erosão típicos conhecidos como piping. 
Ranhuras de escoamento são criadas quando a água sofre erosão e remove material do 
topo do reservatório. Esses sedimentos são transportados para as laterais e para a parte 
inferior, formando pequenos leques. O material retransportado para as laterais forma uma 
nova cobertura de material pelítico sobre o depósito original. As pelotas do depósito original 
diferem daquelas do depósito secundário, pois as primeiras apresentam laminação, enquanto 
as últimas possuem estrutura sólida. 

PALAVRAS-CHAVE: Sedimentos lacustres, Hidroclasticismo, Sulcos de erosão. 


11 INTRODUCCIÓN 


La zona de estudio corresponde a los depósitos lacustres de la Formación El Paso 
(Pleistoceno) ubicados en el piedemonte occidental de las Cumbres Calchaquíes, entre la 
localidad de Cafayate y el paraje Lorohuasi. 

Cafayate se ubica al sudeste de la provincia de Salta, a 1.683 m s.n.m., en las 
coordenadas geográficas 26 04'Latitud Sur y 65º58'Longitud Oeste. Lorohuasi se 
encuentra a 3 kilómetros al sur de la ciudad de Cafayate. Posee una cota de 1.618 ms. n. 
m., y se sitúa en las coordenadas geográficas 26º05'Latitud Sur y 65º58'Longitud Oeste. 
Los depósitos lacustres se encuentran entre los siguientes valores de coordenadas: 26º 04 
Latitud Sur; 65º 52" Longitud Oeste y 26º 06” Latitud Sur; 65º 52" Longitud Oeste. 

Uno de los caminos más rápidos para acceder a la zona de interés, desde la ciudad 
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de Salta, consiste en dirigirse hacia el oeste por la Avenida Juan Domingo Perón/ Ruta 
Provincial Nº28, luego continuar hacia el sur por la Avenida Circunvalación Oeste hasta 
tomar la Ruta Provincial Nº24 en dirección a la localidad de Cerrillos. Desde Cerrillos se 
deben recorrer aproximadamente 150 kilómetros hacia el sur por la Ruta Nacional Nº68. 
Antes de llegar al puente del río Santa María se debe girar hacia el sur y tomar el “Camino 
del Norte”. Este camino, con rumbo norte - sur, es empleado para recorrer los vifedos 
ubicados en el piedemonte occidental de las Cumbres Calchaquíes entre las localidades 
de Cafayate y Tolombón. Por el “Camino del Norte” se deben recorrer 3 kilómetros hacia el 
sur hasta llegar a los depósitos lacustres (Fig. 1). 

En este trabajo se efectuará un análisis geomorfológico para determinar el 
comportamiento de estos sedimentos lacustres y, de las demás unidades geomorfológicas 
identificadas en la zona, frente al accionar de los procesos de meteorización y los agentes 
geológicos. 


3500000 3514000, 


3500000 3507000 


Figura 1. Mapa de ubicación. 
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21 GEOLOGÍA 


Sobre la base de observaciones de campo, el análisis de imágenes satelitales, 
la recopilación bibliográfica y los datos obtenidos de la Hoja Geológica 10e-Cafayate 
(Galván 1981), se determinó que en la zona de estudio afloran las siguiente unidades 
litoestratigráficas: Granito Cerro Amarillo (Carbonífero), Formación Lumbrera (Paleógeno- 
Eoceno), Formación El Paso (Pleistoceno) y depósitos cuaternarios (Fig. 2). 
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Figura 2. Mapa Geológico del área de estudio. 


Carbonífero - Granito Cerro Amarillo (Rapela 1976) 

La roca aflorante corresponde a un granito biotítico-hornbléndico de grano grueso, 
color rosado claro a rosado violáceo en fractura fresca. La característica distintiva es el 
color ocre-amarillento que le otorga la pátina de alteración. En la zona de estudio aflora en 
las laderas occidentales de las Cumbres Calchaquíes (Fig. 2). 


Eoceno. Formación Lumbrera (Moreno 1970) 


Es la unidad superior del Subgrupo Santa Bárbara, Grupo Salta (Cretácico- 
Paleógeno). Está formada por arcilitas y limolitas rojizas, generalmente macizas o con 


bioturbaciones y desarrollo de paleosuelos. Con intercalaciones de arcilitas verdes y grises, 
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calizas amarillentas y areniscas finas. Las arcilitas grises suelen contener altos porcentajes 
de materia orgánica. En el área de trabajo aflora al noreste de las Cumbres Calchaquíes y 
los depósitos corresponden a limolitas y areniscas finas de color rojo y castafio rojizo. 


Pleistoceno. Formación El Paso (Vilela y García 1978) 


Esta unidad está constituida por limolitas y arcilitas pardo claro a pardo amarillento, 
con estratificación fina a muy fina (Fig. 3), intercaladas con areniscas de grano fino. En 
la base son frecuentes las intercalaciones lentiformes de areniscas de grano mediano y 
conglomerados de grano fino. Las pelitas laminadas corresponden a sedimentos lacustres 
que resultan de endicamientos locales relacionados con procesos de remoción en masa 
que, en el Pleistoceno, provocaron interrupciones en los depósitos de abanicos aluviales 
del flanco occidental de las Cumbres Calchaquíes. Las concreciones carbonáticas, rosetas, 
venas y lâminas de yeso fibroso implican condiciones temporales de mayor evaporación en 
el cuerpo de agua. 

En la zona de estudio aflora en el piedemonte occidental de las Cumbres Calchaquíes 
(Fig. 4), formando 7 cordones de relieves de lomadas que varían entre 1,5 a 3 metros de 
altura y 1,5 a 2 metros de ancho. 


Figura 3. Pelitas estratificadas de la Formación El Paso. 


Figura 4. Relieves de lomadas de la Formación El Paso. 
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Holoceno inferior-medio. Cenoglomerados y fanglomerados de abanicos aluviales 


Se denominan cenoglomerados a los materiales depositados por procesos de 
remoción en masa. Sus principales características son: ausencia de estratificación, 
clastos angulares a subangulares, sin selección y matriz fina portante. Los fanglomerados 
corresponden a sedimentos de diversos tamahios, clastos redondeados a subredondeados, 
depositados en un abanico aluvial, por corrientes de alta fluidez y competencia o capacidad 
de transporte (Fig. 5). 

Los fanglomerados y cenoglomerados afloran en el piedemonte occidental de las 
Cumbres Calchaquíes. Somonte y Collantes (2007) a partir del estudio del barniz de las 
rocas, afirman que en el periodo Holoceno existieron ocho eventos húmedos datados entre 
7.300 AP y 300 AP. Estos eventos permitieron la sedimentación de los cenoglomerados y 
fanglomerados. 


Figura 5. Fanglomerados intercalados con sedimentos más finos. 


Holoceno inferior-medio. Gravas, arenas y limos fluviales 


Los depósitos fluviales en la zona corresponden a limos, arenas y gravas que forman 
parte de las terrazas fluviales, la Ilanura de inundación y el cauce actual del río Santa María 
(Fig. 2). 

Las terrazas fluviales son discontinuas a lo largo del río Santa María y en algunos 
sectores presentan una capa de arena rica en materia orgânica datada en 2.190 + 530 AP 
por Strecker (1987, citado por Collantes y Busnelli, 2014). 


Holoceno inferior-medio. Arenas y limos eólicos 


Desde Cafayate hasta Lorohuasi, se extienden extensos campos dunarios, actuales 
y pasados, sobre parte de la Ilanura de inundación y terrazas fluviales del río Santa María 
(Fig. 6) y, sobre los extremos medios-distales de los abanicos aluviales ubicados en el 
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piedemonte occidental de las Cumbres Calchaquíes (Fig. 2). 

Estas dunas están constituídas por arenas de color gris pardo, con muy pobre 
selección eólica, granos muy finos a medianos, en su mayoría, subredondeados a 
subangulosos. La matriz es limosa y de color castafio claro. La mineralogia corresponde 
a cuarzo, feldespato potásico y plagioclasa, con pequefias cantidades de biotita, clorita 
y cristales volcánicos. Además, las arenas contienen una gran variedad de minerales 
pesados, con augita como principal componente, junto con zircón, turmalina y horblenda 
(Pefia Monné et al., 2016). 

De acuerdo al estudio realizado por Pefia Monné et al. (2016), durante el último 
milenio, las fases húmedas y secas alternantes favorecieron la movilización y estabilización 
de las dunas. Estas fases se establecieron mediante dataciones OSL (dataciones por 
Luminiscencia Ópticamente Estimulada), restos arqueológicos y datos históricos. Los 
intervalos secos ocurrieron entre ca. 1.000-1.100 AP, 1.300-1.420 AP, 1.550-1.680 AP, 
1.740-1.850 AP y el período cálido actual. 


DD "aa 


Figura 6. Terraza parcialmente cubierta por depósitos eólicos actuales, margen oeste del rio Santa 
Maria. 


31 GEOMORFOLOGIA 


Una Unidad Geomorfológica es una forma característica del relieve, definida por 
una génesis particular, además de ser cartografiable en un intervalo de escala establecido 
(Sayago, 1982). Siguiendo los criterios establecidos por Van Zuidam y Van Zuidam 
Cancelado (1979), Verstappen y Van Zuidam (1991) y Collantes (2001), en la zona de 
estudio se reconocieron desde el punto de vista morfogenético: 1) Unidades de Origen 
Estructural-Denudativo, 2) Unidades de Origen Fluvial-Aluvial, 3) Unidades de Origen 
Fluvial y, 4) Unidades de Origen Eólico (Fig. 7). 
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Figura 7. Mapa Geomorfológico. 


Unidades de Origen Estructural-Denudativo 


Se agrupan en esta clase las formas desarrolladas bajo la influencia de factores 
geológicos endógenos (litologíia y estructuras) y de procesos exógenos (erosión, 
meteorización, remoción en masa, entre otros). En el área de estudio, corresponden a esta 
categoría, las laderas occidentales de las Cumbres Calchaquíes y los relieves de lomadas 
de la Formación El Paso. 

Las laderas constituyen una importante Unidad Geomorfológica desde el punto de 
vista espacial, caracterizadas en general por un relieve abrupto, con pendientes mayores a 
los 20º, con desarrollo de valles profundos y empinados, escarpas estructurales y cicatrices 
de deslizamientos actuales y fósiles. Están afectadas por: 1) alta densidad del drenaje, 2) 
meteorización física y, 3) procesos de remoción en masa actuales y pasados, que en la 
zona adquieren magnitud debido a la pendiente, las precipitaciones pluviales torrenciales 
y la escasa vegetación. 

Los principales procesos de meteorización física que actúan sobre las laderas están 
determinados por las variaciones térmicas propias del clima árido de la zona y, corresponden 
al termoclastismo y crioclástismo. 

Las laderas están expuestas a la insolación diaria, Ilegando a alcanzar una 
temperatura media de 30ºC, esto combinado con el enfriamiento nocturno, hace que la 
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variación de temperatura diaria varie entre 6º y 30ºC, por lo que las rocas experimentan 
dilatación y contracción, favoreciendo eltermoclastismo. El crioclástismo actúa, en la cima de 
los relieves positivos, durante la estación invernal. Se desarrolla a partir del congelamiento 
y descongelamiento del agua en las diaclasas y fracturas de las rocas. Como consecuencia 
de ambos procesos se origina una importante cantidad de material detrítico. 

Por otro lado, en el piedemonte occidental de las Cumbres Calchaquíes, se 
reconocieron relieves de lomadas correspondientes a afloramientos de los depósitos 
lacustres de la Formación El Paso. Estos relieves son el resultado de la intensa erosión e 
incisión fluvial, que afecto la continuidad del primitivo deposito lacustre, del cual se pueden 
ver afloramientos remanentes, entre los abanicos aluviales. 

En los depósitos lacustres el escurrimiento concentrado genera piping (Fig. 9). Estas 
morfologías son exclusivas de las regiones secas. Los surcos de escurrimiento se originan 
cuando el agua erosiona y elimina el material que se encuentra en la parte superior del 
depósito. Este material es transportado hacia los laterales. De esta manera se forma sobre 
el depósito original una nueva cobertera de material pelítico. Las pelitas del depósito original 
se diferencian de las del depósito secundario, porque las primeras presentan laminación, 
mientras que las segundas tienen una estructura maciza. 


Figura 9. Piping en depósitos de la Formación El Paso. 


Unidades de Origen Fluvial-Aluvial 


De acuerdo a Verstappen (1983), en este grupo se incluyen todas aquellas formas 
en cuya génesis intervienen la dinâmica fluvial y los procesos de remoción en masa. 

En el área de trabajo estas Unidades se encuentran en el piedemonte occidental de 
las Cumbres Calchaquíes y están representadas por abanicos aluviales (Fig. 7). 

Desde el punto de vista de su desarrollo morfogenético, se identificaron tres ciclos 
principales de aluvionamiento, dos que constituyen la porción inactiva y uno correspondiente 
al área activa (Fig. 7). El primer ciclo, ubicado en el extremo septentrional, es el más 
antiguo, corresponde a remanentes de la depositación que habría dado nacimiento a los 
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abanicos, presenta como rasgo sobresaliente la presencia de paleocauces que recorren 
toda su superficie. El segundo ciclo esta caracterizado morfológicamente por abanicos más 
bien restringidos, afectados intensamente por la erosión hídrica (formación de surcos y 
cárcavas), desarrollados sobre parte de la superficie correspondiente a los abanicos del 
primer ciclo. El tercer ciclo está constituido por los abanicos actualmente activos, dichos 
abanicos se caracterizan por ser de gran extensión, unirse unos con otros en sus extremos 
distales y ser, de los tres ciclos mencionados, los menos afectados por los procesos de 
erosión hídrica y los más afectados por la actividad eólica. 


Unidades de Origen Fluvial 


Se consideran Unidades de Origen Fluvial aquellas geoformas en cuya génesis la 
dinámica fluvial es dominante. En el área de estudio esta clase está representada por el 
valle fluvial del río Santa María: cauce actual- Ilanura de inundación y terrazas fluviales. 

En la zona de estudio el rio Santa María constituye el nivel de base local. El cauce 
actual presenta un disefio de tipo fluvial entrelazado. Este diserio se origina debido a que el 
curso fluvial no puede transportar toda la carga que recibe y en consecuencia la abandona 
en su lecho, formando barras e islas. Estas barras van dividiendo la corriente en cauces 
secundarios y en el verano, en períodos de crecidas, pueden quedar sumergidas por el 
agua. El exceso de carga que recibe el rio Santa María proviene: 1) de ríos tributarios, los 
cuales durante el verano se ven afectados por crecidas al ocasionarse Iluvias torrenciales y 
2) de los procesos de remoción en masa, principalmente flujos densos, el material detrítico 
transportado por los procesos gravitacionales proviene, fundamentalmente, del accionar de 
la meteorización física sobre las rocas de las laderas de las Cumbres Calchaquíes. 


Unidades de Origen Eólico 


En el área de estudio el clima árido favorece el accionar del viento como un agente 
geomorfológico que erosiona, transporta y deposita. Se reconocieron, en base a las 
observaciones realizadas en el campo, el análisis de imágenes satelitales y la recopilación 
bibliográfica: dunas linguoides-barjanoides, nebkas y dunas longitudinales-lineales. 

Las dunas linguoides—barjanoides son dunas complejas, con una forma intermedia 
entre barjanes aislados (dunas individuales con formas de media luna, con sus extremos 
apuntando en la dirección del viento) y dunas transversales (dunas compuestas, constituídas 
por cordones fundamentalmente perpendiculares separados por depresiones orientados 
con ângulos rectos respecto al viento que posee dirección e intensidad constantes). 

Las nebkas son dunas simples, obstaculizadas por la vegetación. No superan el 
metro de altura y se forman a la sombra de un obstáculo, como una mata vegetal o una 
roca. 

Las dunas longitudinales o lineales son dunas complejas, representadas por crestas 


rectilíneas, paralelas, con un espaciado regular y una sección transversal más simétrica 
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que las dunas transversales. Se forman y se desplazan paralelas a la dirección del viento 
predominante a la que se superponen vientos oblicuos que arrojan arena lateralmente. 
Requieren de un limitado abastecimiento de arena y de un terreno duro y rocoso. 


41 CONCLUSIONES 


Teniendo en cuenta el clima árido del tipo BWk - según la clasificación de Kôppen- 
Geiger (1918 y 1936) en Kottek et al. (2006) - y el análisis geomorfológico efectuado, en 
la zona de estudio, el principal agente geomorfológico es el agua y está relacionado con 
las precipitaciones medias anuales; para la zona se estima un valor medio de 204 mm/ 
anuales, de acuerdo a los datos proporcionados por la Bodega La Rosa, para el período 
2000-2010. Además, las precipitaciones se caracterizan por su torrencialidad, lo que da 
lugar a una intensa erosión hídrica pluvial y fluvial. El otro agente es el viento. 

El accionar de la erosión hídrica se manifiesta en el desarrollo de surcos y cárcavas 
sobre los abanicos aluviales del piedemonte occidental de las Cumbres Calchaquíes y en 
el desarrollo de piping en los depósitos lacustres de la Formación El Paso. 

El viento actúa erosionado los arbustos de la zona, vifiedos, algarrobos, entre 
otros, dejando expuestas sus raíces y, transportando y depositando arena, desarrollando 
extensos campos dunarios. 
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